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ABSTRACT 



The fundamental concept of orthogonality of mathematical objects occurs in a wide 
variety of physical and engineering disciplines. The theory of orthogonal functions, for 
example, is central to the development of Fourier series and wavelets, essential for signal 
processing. In particular, various families of classical orthogonal polynomials have tra- 
ditionally been applied to fields such as electrostatics, numerical analysis, and many 
others. 

This thesis develops the main ideas necessary for understanding the classical theory 
of orthogonal polynomials. Special emphasis is given to the Jacobi polynomials and to 
certain important subclasses and generalizations, some recently discovered. Using the 
theory of hypergeometric power series and their q —extensions, various structural prop- 
erties and relations between these classes are systematically investigated. Recently, these 
classes have found significant applications in coding theory and the study of angular 
momentum, and hold much promise for future applications. 
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I. INTRODUCTION 



The abstract concept of orthogonality of functions (or other mathematical objects) 
is a generalization of the notion of having two or more vectors perpendicular to one 
another. This concept arises naturally in a wide variety of physical and engineering 



ment of Fourier series and wavelets which are essential to signal processing. 

Classical Fourier series (real form) depend on the property that the trigonometric 
functions sine and cosine are orthogonal (on an appropriate real interval) in a formal 
sense that will be made precise later. As a consequence, a bounded periodic function 
f(x) of period 2 n which satisfies the Dirichlet conditions! may be expressed in the form 



are the classical Fourier coefficients. These formulas can be modified via a change of 
variable to accomodate any such function of period 2 L. (Ref 1: p. 529] 

This property can be used to generate other classes of orthogonal functions - 
polynomials, for example - that behave in very structured and useful ways such as in 
generalized Fourier series. In particular, specific families of these "classical orthogonal 
polynomials" have traditionally been used for solving problems arising in various areas 
of applied mathematics, physics, and engineering, among others. 

This thesis develops the main ideas necessary for understanding the classical theory 
of orthogonal polynomials. Special emphasis is given to the Jacobi polynomials and to 
certain important subclasses and generalizations. Much of the investigation will be 

1 Dirchlet conditions: (i) In any period f(.x) is continuous, except possibly for a finite number 
of jump discontinuities, (ii) In any period f(x) has only a finite number of maxima and minima. 



disciplines. For example, the theory of orthogonal functions is central to the develop- 



OO 




where 
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made using the theory of hypcrgeomctric power series and their q —extensions. The 
classes discussed in this thesis arc but a small fraction of those identified and studied in 
the literature. 

A. CHEBYSHEV POLYNOMIALS 

The Chebyshev polynomials of the first kind, T n {x), arise from an elementary trigono- 
metric consideration. As such, they satisfy various properties and identities which arc 
easily derived directly from their definition, many of which arc observable from their 
graphs (see below). This class of polynomials will serve as the model for some of the 
basic structure of more general classes. 

For n = 0,1,2, ... , define 

T n {x ) = cos(/7 arccos x), — 1 <; x < 1 

i.e., letting x = cos 0 , 0 <; 0 <, n, 

(1) T n { cos 0) = cos n0. 

Some immediate consequences of (1) are |r„(ar)|<l for |x| < 1, with 
r n (cos- 7 p-) = (-l)*, 0 < k < n ; in particular T n (l)= 1 and T n ( -1) = ( -1)" for all n 
which can be seen graphically in Figure 1. 

1. Three-term Recurrence Relation/Differential Equation 

From (1), we have 



(2) T 0 (x) = 1 and r,(x) = x, 
and by considering the identity 

(3) cos(a + b) + cos(a — b) = 2 cos a cos b 
with a = nO, b — 0, we obtain 

W T n+l (x) = 2xT n (x) - T„_ t (x). 

Equation (4) is known as the three-term recurrence relation for T^x) which to- 
gether with initial conditions (2) imply that T„{x) is a polynomial of degree exactly n, 
called the n"‘ Chebyshev polynomial of the first kind. Note that the leading coellieicnt 
of T n (x) is 2' ,_1 for n > 2. An inductive argument applied to this recursion shows that 



2 




Figure 1. Chebyshev Polynomials 

T„(x) is an even function if n is even, and odd if n is odd (see Figure 1). The first few are 
listed below: 



T 0 (x) = 1 , T l (x) = x, T 2 (x) = 2* 2 - 1 , T 3 (x) = 4r 3 — 3.v 

T 4 (x) = 8x 4 — 8 jc 2 + 1 , T 5 (x) = l6x S — 20x 3 + 5x. 

Differentiating (1) twice with respect to 0 yields the second order differential 
equation for T n {x) : 

(5) (i - * j ) 7VM - *7VM + * J r„w = o. 
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2. Orthogonality of Chebyshev Polynomials 
Let 



Jo 



I mn = I cos mO cos nOdO. 



If m # n, then using (3) with a = m0, b = nQ yields 

1 f sin(m + ri)Q sin(/n — n)Q 



=-T 

n ' n 2 l 



m + n 



+ 



m 



^1 - 



0. 



If m = n 0, then by using the identity cos 2 a = — (1 + cos 2a), 



/ =-r 

n ' n 2 L 



d + -r— sin 2n0 

2n 



7T 

] -f- 

— 2 n 



If am = n = 0, then 



3,0 = [ dO = 7Z. 



Hence, 



( 6 ) 

where 



and 



f 



cos cos dO = h n x 5 mn 




n, m = n = 0 
7r/2, m = n ^ 0 



<3 






0, rn # n 

1 , /n = n 



is the Kroncckcr delta function. Changing variables via jv = cos 0, we have 
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(7) 



J 



TJx)TJx) 



■ dx = h„ 1 3 



v“ 



mjv 



This important property is formally known as the orthogonality relation for the 
Chcbyshev polynomials. The reason for this terminology will become clear in the next 
chapter. 

3. Zeros of Cliebyshev Polynomials 

2k — 1 

Setting T n ( cos 9) = cos nO = 0, we obtain 9 = 9 k n = — - n, i.c., 

In 

X = X*, „ = cos 0 K „ , 1 < k < n. 

Thus all the zeros of T n {x) are real, distinct, and may be regarded as the 
projections onto the interval (—1,1) of the equally distributed points 9 k n on the unit 
circle, as seen in Figure 2 . Moreover, for 1 < k < n, an easy algebraic check verifies that 
Ok, „+■ < 0*, n < 0 M . *+i and therefore x K n+1 < x k<n < x* +li n+1 . Hence, the zeros of T n+l (x) 
interlace with those of T n (x). This interlacing of zeros is a striking feature of the plots 
in Figure 1. 

The zeros of Cliebyshev polynomials, and of other orthogonal polynomials in 
general, are extremely important for applications to numerical analysis, electrostatics, 
and many other fields. 

4. Looking Ahead 

Many of the properties derived for the Cliebyshev polynomials T„(x) from their 
trigonometric definition (1), extend to more general classes of orthogonal polynomials 
via a general theory, elements of which will be developed in this thesis. Some of the 
many properties satisfied by these classes that we will derive include: 

1. Orthogonality with respect to a weight function 

2. Three-term recurrence relation 

3. Second order differential or difference equation 

4. Hypergeometric series expression 

5. Rodrigues' formula 

6. Generating function. 

The general approach we will take is to define these "classical orthogonal polynomials" 
via terminating hypergeometric power series, and from this prove (most of) the other 
properties. However, because of this equivalence, many authors choose to define a given 
class using one of these other characterizing properties. 
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Figure 2. Zeros of Chebyshev Polynomials 



In order to understand the general theory, it is first necessary to define the ab- 
stract concept of orthogonality in an appropriately defined "space" of functions. We 
turn our attention to these fundamental ideas in the next chapter. 
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II. BACKGROUND 



A. ELEMENTARY LINEAR ALGEBRA 
1. Vector Spaces 

Let R" denote the collection of all vectors ( 11 -tuples), u = (<7,, a 2 , ... , a n ), where 
each a, e R, /= 1, 2, ... , n. The standard Euclidean inner product (also referred to as the 
dot product ) of two such vectors u = (a„ a 2 , ... , a n ) and v = {b u b 2 , ... , b„ ) is given by 



Two vectors u, v e R" arc perpendicular, or orthogonal, if and only if < u, v> = 0. 

The objective of this chapter is to extend these familiar notions to objects other 
than classical Euclidean vectors, in particular, the "vector space" of polynomials defined 
on a real interval [a, b~\. 

A vector space V over a scalar field F (usually R or C ) is a nonempty set of 
objects called vectors, for which the operations of addition and scalar multiplication arc 
defined. Addition is a rule for associating with each pair of vectors u and v in V an 
element u + v, called the sum of u and v . Scalar multiplication is a rule for associating 
with each scalar c in F and each vector u in V an clement cu, called the scalar multiple 
of u by c. [Ref. 2: p. 150] 

For all u, v,we V and c,d e F, a vector space V must satisfy: 

1. Additive closure. 11 , v e V => u + v e V 

2. Commutativity, u + v = v + u 

3. Associativity. 11 + (v + w) = (u + v) + w 

4. Additive identity. There exists a zero vector, 0 e V, such that 0 + u = u 4- 0 = u. 

5. Additive inverse. For each u e V, there exists a vector — u e V, such that 
u + ( - u) = ( - u) + u = 0. 

6. Multiplicative closure, ue f and c e /•' => cu e V 



n 




The length, or norm, of a vector u e R" is given by 
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7. Distributivity. c(u + v) = cu + cv 

8. Distributivity. (c + d ) u = cu + du 

9. Multiplicative associativity. c(c/u) = (cd ) u 

10. Multiplicative identity. There exists a scalar 1 e F such that lu = u. 

Example 1: FT (the model) 

Example 2: Pf_a,d] = (polynomials of degree < N on the interval [a,/;]} 

Example 3: P\_a,d] = (polynomials on the interval [a, d]} 

Example 4: C\_a,b ] = (continuous functions on the interval [a,/;]} 

Note that P N [.a,d] c P\_a,d] c C\_a,b~]. These vector spaces are sometimes re- 
ferred to as function spaces. The interval [a, 6] may be finite or infinite (i. e., 
\_a, oo) , (— oo, d] , or (— oo, oo) ) for our purposes. 

A subset U of a vector space V is said to be a vector subspace of V if it is a 
vector space in its own right. 

Example 5: Pf_af] is a subspace of P\_a,d], which in turn is a subspace of 

Cla,b}. 

Given a set of vectors (v„ v 2 , ... , v„] in a vector space V, and scalars 
c„ c 2 , ... , c„, the vector c,v, + c 2 v 2 + ••• + c„v„ e V is said to be a linear combination of 
(v,, v 2 , ... , v„}. The set of vectors (v,, v 2 , ... , v„} is said to be linearly dependent if there 
exist scalars c„ c 2 , ... , c„, not all equal to zero, such that the linear combination 
c,v, 4- c 2 v 2 + ••• + c n \ n = 0. (Equivalently, at least one of the vectors v, can be expressed 
as a linear combination of the others.) Otherwise, the set (v,, v 2 , ... , vj is linearly inde- 
pendent. An infinite set 5 = (v„ v 2 , ... , v n , ... } is defined to be linearly independent if every 
finite subset of 5 is linearly independent; otherwise S is linearly dependent [Ref. 3: p. 8], 
The vectors (v t , v 2 , ... , v„) are said to span V if every vector \ e V can be represented as 
a linear combination of (v„ v 2 , ... , v„}. In this case, we write V = span{v t , v 2 , ... , v„}. 'fhe 
vectors (v,, v 2 , ... , v„} form a basis for V if they are linearly independent and span V. The 
dimension of V is the number of elements in any basis. 

Example 6: The set (e,, e 2 , ... , e„} is the standard basis for R", where 

e, = (0, 0, ... ,0, 1,0, ... , 0) i.e., the vector with a one in the i' h position and zeros else- 
where, / = 1,2, ... , n. 

Example 7: The set (1, x, .x 2 , ... , .v ,v } is the standard basis for P Ja.b']. (Linear 
independence is ensured by the Fundamental Theorem of Algebra.) The dimension of 
Pf_a,b~] is therefore iV+1. 
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Example 8: The set {1, .r, x 2 , ... , x", ... } is the standard basis for P\_a,b~\, and 
hence P[_a,b'] is an infinite-dimensional vector space. 

2. Inner Product Spaces 

An inner product on a real vector space V is a mapping 



such that for all u, v, w e V and a,P eU, the following properties hold: 

1. Positive definiteness: < u, u> > 0, and < u, u> = 0 if and only if u = 0 

2. Symmetry: < u, v> = <v, u> 

3. Bilinearity: < au + fix, w> = a < u, w> + /? < v, vv> 

A vector space with an inner product is known as an inner product space. 

Example 9: V = IR" ; let constant "weights" w ( > 0 be given, / = 1,2, ... , n. 
For u = (a,, a 2 , ... , a n ) and v = (&„ b 2 , ... , b n ), u, vef, 



If w,= 1 for /= 1,2,... , n, then this reduces to the standard Euclidean inner product, 
or dot product. Otherwise, this is referred to as a weighted inner product. 

The next two examples are commonly applied inner products on function space, 
and are analogues of the previous example. We assume a given weight function 
vv(a-) > 0 in {a,b), integrable in the first case (e.g., continuous for [afi] a finite interval). 



< , >:VxV^U 



n 




Example 10: V = P N [a,b~\, P\a,b], C\_a,b~\ 



<fi,g> = f(x)g(x)w(x)dx 






Example 11: V = Pfiafi] 




(Positive definiteness is ensured by the Fundamental Theorem of Algebra.) 



The norm induced by the inner product is given by ||u|| = 7 < n,u> .2 
Example 12: For the inner products of Examples 10 and 1 1 therefore 



ll/ll = 




and 



( N 

ll/ll = ( X! [/W ] 2 “'W 

V x=0 

respectively. These are sometimes referred to as "Z, 2 -norms." 

Two vectors u,ve f are said to be orthogonal, denoted u_Lv, if and only if 
<u, v> =0. The vectors u and v are said to be orthonormal if ulv and 
INI = ||v|| = 1. Note that the orthogonality of vectors in a space is determined by the in- 
ner product being used. 

The two examples which follow refer back to Chapter I, Section A.l. 

Example 13: Formula (6) shows that the functions { 1, cosx, cos 2x, ... } are 
orthogonal on [0, 7r] with respect to the uniform weight function vv(x) = 1. A similar 
computation shows that the same property holds on [ — tt, n ] with respect to the weight 
function tv(.v) = -jp , i.e., 




</,£>= 4* f fix) Six) dx. 

One advantage of preferring this inner product over the standard one lies in the com- 
putation of norms. Using ||/|| = J < f, f > , we have || 1 || = 2 and || cos nx || = 1 if 
n> 1. Hence the functions { 1/2, cosx, cos 2x, ... } are orthonormal on [— n, rc] with 
respect to the inner product above. Similar statements hold for the integral of a product 
of two sine functions on [—n, re], as well as for the product of a sine and a cosine. 

Example 14: By (7) the Chcbyshev polynomials {7’„(x)} form an orthogonal 

class with respect to the inner product of Example 10 above on [ —1,1] with the weight 
function w(x) = (1 — x 2 )" 12 . 

2 Recall that < u, u > >0. We remark that in the same way we defined inner product earlier, 
it is possible to define a general norm on a vector space wliich is not induced by an inner product. 
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Wc remark here that Examples 10 and 11 can be unified into a single inner 
product on a "polynomial" space V via 



where da(x) is a positive Lebesgue-Stieltjcs measure on a measurable set E possessing 
finite moments, i.e., x n da.(x) integrable, n = 0, 1,2,.... In Example 10, E = \_a,b] <= R and 
da.(x) = vv(.r)^.v; the resulting expression is known as a continuous inner product, while 
in Example 11 the set £ consists of a finite number of points {0, 1, ... , N] cR, and the 
associated measure gives rise to a discrete inner product. 

B. FOURIER SERIES 

Let v e V, and U be an Aj-dimensional subspace of V having some orthonormal basis 
{ u,, ... , u„ }. (Any basis can be orthonormalized via the Gram-Schmidt process - see next 
section.) The vector v can be resolved into a sum of two vector components: 

(1) y = (v — w) + w 

where w e U and (v — w) A. U. (See Figure 3.) The vector w is referred to as the 
orthogonal projection of v onto U. Since the vector (v — w) is orthogonal to every vector 
in U by construction, it follows that for each j — 1,2, ... , n, < v— w , u, > = 0, or 

(2) < V, Uy > = < W, Uy>. 



Moreover, since it lies in U, vector w can be expressed as some linear combination of 
{ u„ ... , u„ } : 



Take the inner product of both sides with u y for each j — 1,2, ... , n. From the assumption 
that < u, , u y > = 0 unless i = j, we have the property that 




E 



n 




< W , Uy > = Cj < Uy , Uy > . 



Thus, 



Cj = < v , u ; > 
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via (2) and the assumption that < u, , u, > = || u ; || J = 1. 

Thus, 

n 

(4) w = Yj < v > u < > u ' 

1=1 

and this vector represents the "best approximation" in U to v e V in the sense that of 
all vectors z e U, it is the projection vector w e U which uniquely minimizes the distance 
II v-z || . 

Suppose now that U is an infinite-dimensional subspace of V (also infinite- 
dimensional), having orthonormal basis { u„ ... , u„, ... }. Then from (1) and (4) 

n 

V = (v-w) + < V , u ; > u f 

1=1 



we may write 



(5) v = X<v,m>u, 

;=i 

in the sense that 

(6) lim || v— w || = 0, 

n—*oo 

i.e., the norm of the "residual vector" (and hence the vector itself) v— \v -» 0 as 
— > co. Formula (5) is known as the generalized Fourier series for v e V with respect 
to the orthonormal basis { u, The coefficients given in (3) arc called the generalized 
Fourier coefficients of v e V. Statement (6) is known as the norm convergence property 
of Fourier series, and the "minimization property" mentioned above extends to this 
infinite-dimensional case. 

Example 15: Let V= C[a,6], and { 0,(.v) }" 0 be an orthonormal basis of 

eigenfunctions (sometimes referred to simply as an eigenbasis ) of V. Then f e V has a 
Fourier scries representation 



f(f ~ c ‘ ( W- V ) 

7=0 
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V 




with Fourier coefficients 



(7) c t = < f, <f>i > = f f(x) <p t (x) w(x) dx. 

^ a 

In this function space context, norm convergence 



lim 

n-+oo 




n 



fix) - Yj c l 
/= 0 
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is referred to as mean square convergence, and is the least squares principle in regression 
analysis. 

In particular, if V is equipped with the inner product of Example 13 and 
orthonormal basis 



on [—7i, 7r] (see Example 13, Section A. 2), then a suitable function feV (and its 
2/T-periodic extension on IR ) has a classical Fourier series 



as indicated in the Introduction. 

C. GRAM-SCHMIDT ORTHONORMALIZATION 

The Gram-Schmidt process orthonormalizes any set of linearly independent vectors 
in an inner product space. This method will be used in later sections for different inner 
products on the vector space P\_a,b]. 

Begin with an inner product space V and any set of vectors {v,, v 2 , ... , v„, ... }, finite 
or infinite, such that any finite number of elements of this set are linearly independent. 
Recursively define a new set of vectors {u„ u 2 , ... , u„, ... } 



{ 4>i(x ) } = { 1/2, cos .r, cos 2.r, ... , sin x, sin 2x , ... } 




where 



1 f * 

a n = < f, cos nx > =— f{x) cos nx dx 



K - </, 





, k 1 , 2 ,...,/?,... 



where y* = v k — n*, with3 



3 By convention, V<j, = 0, giving \v, = 0. 
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These new vectors {u„ u 2 , ... , u„, ... } are orthonormal by construction and span the 
same space as the original vectors. Note that this process occurs in two stages: 
orthogonalization and normalization. The orthogonalization is accomplished by sub- 
tracting \v k , the orthogonal projection of \ k onto the subspacc spanned by 
{u,, u 2 , ... , u*_,}. The component of v* which remains, denoted above as y k , is then 
orthogonal to the vectors {u,, u 2 , ... , u*_,} as shown in Figure 4. The normalization is 
then achieved by dividing y* by its norm, thus giving it unit "length". 

1. Legendre Polynomials 

Example 16: Let V = P\_ — 1,1] with basis {1, .r, x 2 , ... , x *, ... } and uniform 

weight function w(x) = 1. The inner product is then given by < f,g > = f'^f(x)g(x)dx. 
The Gram-Schmidt process yields the set 




as an orthonormal basis for P[ — 1,1]. Since this set is linearly independent, we can 
standardize the set by taking scalar multiples of these polynomials so that P n (l)= 1. 
Members of the resulting orthogonal set 



are known as the Legendre polynomials on [—1,1]. If the normalized Legendre 
polynomials { p„{x) }~ 0 are used as the orthogonal eigenbasis for a Fourier series, the re- 
sulting expansion is often referred to as a Legendre series representation; when 
Chebyshev polynomials are used, we obtain a Chebyshev series representation, etc. 

The Gram-Schmidt process can always be used in this way to generate a class 
of orthogonal polynomials with respect to a given inner product (i.e., weight function) 
on a real interval. When using the Gram-Schmidt process from the basis 
{l,ar, .v 2 , ... , .v", ... }, the orthogonalization stage producing y k results in a set of monic 
polynomials, i.e., the leading coefficient of each polynomial is one. In the normalization 
stage, we are dividing by the norm || y k || > 0. Thus the leading coefficient of polynomials 
in an orthogonal class is strictly ‘positive. In the next chapter, we will examine other 
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spanluj ,112) 




Figure 4. Orthonormalization 

ways to define these classes. It is the structure and applications of certain of these 
classes with which we will primarily be concerned. 

D. THE GAMMA FUNCTION 

The gamma function r(.c) is a fundamental mathematical object that appears fre- 
quently in the representations of orthogonal polynomials as well as in many other ap- 
plications. This "special function" was developed as a generalization of the factorial 
function of the natural numbers. As we will see, the gamma function has the value 
(n — 1)! for the positive integers n but it is defined for noninteger values as well. 

A conventional definition for the gamma function is 
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( 3 ) 



F(-r) = 



e ! t x l dt , x > 0. 

•'o 

The positivity of jc ensures that this improper integral converges. We now develop some 
fundamental properties of the gamma function. Integration by parts in (8) yields 

(9) r(x+ 1) = jcr(jc). 

We now introduce the Pochhamtner symbol or shifted factorial, («)„ , to simplify our 
notation. For n> 0, define 

( a) n = a(a + l)(n + 2) ... (a + n — 1) , if n > 1 

and (a) 0 = 1. Letting a = 1 gives (1)„ = (1)(2)(3) ... («) = n\. Note that for a negative inte- 
ger, (— m) n = 0 if n> m> 0. The shifted factorial can be defined for negative subscripts 
but we will not need this in our work with polynomials. 

Iteration of (9) n times yields 



(10) r(.r + «) = (.v)„rw 

for every positive integer rt. Using this property, the gamma function can be extended 
to include negative real numbers by defining 

(11) r(.r)=— — r(x + /z) for -«<x<-/7+l. 

UN 

Since this expression is undefined when jc is zero or a negative integer, the gamma 
function is not defined for those values. 

Letting x = 1 in (8) and computing directly, we have T(l) = 1. It then follows that 
r(n + 1) = n\ by letting x = 1 in ( 10). Furthermore, 



/ . \ 


p°° 


T = 


e-'r^dt = 2 


\ Z / V 


o ^ 



Ooc 




where the second integral can be evaluated by standard methods involving multiple in- 
tegrals. 
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Finally, we define a generalized binomial coefficient as follows, 
negative integers, define 



For .v and a non- 



/.r+a\ f.r+a)! 

{ x x\a\ ’ 

For nonintegral a, define 

/x+a\ (a+l)jc r(jr+a+l) 

V X ) = (l) x = r(jr+I)r(a+l) • 



1. The Beta Function 

An integral related to the gamma function defines another useful function called 
the beta function which is given by 



(12) B(x,y)= f t x -\\-tY~ x dt 

n> 

for jr,,y>0. We now establish an important connection between the beta and gamma 
functions. We start with an identity easily verified from (12): 

(13) D(x, y + 1) = D(x, y) - B(x + 1, y). 

Also from (12) 



B(x+l,y) = 



r T (i - t) y ~ x dt 



‘'O 



which when integrated by parts gives 



B(x+ 1, y)=yB(x, y+ 1). 
Substituting into (13), we obtain 

B(x, >0 = — B(x, y + 1) 



which when iterated yields 
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, (*+y)n n/ . , (*+A 

D ( x ’ y) = “ 73 — 5 (- v ’ >' + ") = • 



(A 1 ’ ' 7 (A 

Changing variables from t to — , 



t x ~\i _ t f+y-'dt. 



( 14 ) 



D{x, y ) = 



(*+A 
(A n x 



r 



x-1 



(l-T 



Taking the limit as n 



oo , and using the fact that lim(l — — = <r f , we have 



(15) 



B(x, y) = r ( ,)l i m 



C*+A 
(A«* ' 



(The fact that we can pass the limit through the integral on the right can be mathemat- 
ically justified.) I fy = 1 , then (15) gives 



B{x,l) = T(jc) lim 

n-*oo 



(•* + i)^i 
n\ n x 



By direct evaluation using (12), 



5(*,1) 



f 1 

4 



Hence 



= T(.t) lim 

n — * oo 



(•* + On 

I * 

n\ n 



which can be written as 



T(.v) = lim 

n-* oo 



n\ n x 

x(x+ 1)„ ' 



Noting that .v(.v + 1)„ = (.t) n+ , = (.v)„(.v + //), we have 
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. l 
n\ n n 



r(jc) = Jim 

n~> oo 



(AC* + n ) 



which gives the form 

(16) r» = lim 44- 

(*)„ 



(Equation (16) was Euler's original definition of the gamma function. A separate "esti- 
mation" argument may be used to show that this limit mathematically exists.) Thus by 
(15), 



D(x, y) = I» lim 

n—+ oo 



(*+A 

7! 

(y) n 



Then by (16), we have the useful identity 



(17) 



B(x, y) 



rwr(y) 

r(x+y) • 



We will find this identity useful in understanding the Jacobi polynomials in 
Chapter IV, where it becomes necessary to evaluate a related integral: 



f ( 1 — xY ( 1 -hr/ dx. 

J -i 

We remark here for future reference that the formal change of variable .r = 1 — 2/ can be 
used to transform this integral into 



2 * +p+ 1 [ 1 /“ (\—tf dt = 2 ,+ ^ +1 B{a+ 1, p+ 1) 

_ ~*+/?+i rfo+O r(/?+i) 

T(a+/?+2) 
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III. GENERAL THEORY OF CLASSICAL ORTHOGONAL 

POLYNOMIALS 

In this chapter, we examine some of the characteristic properties associated with 
classes of orthogonal polynomials. Some of these properties provide alternate means 
of defining a class. These alternate definitions often provide a straightforward way of 
producing a specific result that may be very difficult to derive otherwise. 

Throughout this chapter, we let { />„(. r)}~ 0 denote a set of real polynomials with 
p n (x) of degree n, i.e., 

(1) p n {x) = k n x n + s n x n ~ 1 +••• , /<„ > 0. 

Recall that these polynomials are said to be orthogonal on an interval [a, 6] with respect 
to a continuous weight function w(jc) > 0 on {a, b) if 

(2) < Pm - Pn > = [ Pm( x )Pn( x ) w ( x ) dx = h~ X 5 m n 

where the normalization h~ l # 0 is chosen to simplify the expression of certain formulas. 
Note that since 

II Pn II 2 = f [ Pn( x ) fw{x)dx = h ~ x , 



it follows that h„ > 0. 

A. POLYNOMIAL EXPANSIONS 

We begin by showing that any real polynomial ^ m (.xr) of degree m on [a,/;] can be 
written as a linear combination of orthogonal polynomials { /? n (.v)}” =0 : 



m 

(3) q m (x) = Yj °Pm Pk x ) 

i=0 

for constants ot, m , / = 0, 1,2, ... , m. 

The proof is by induction on the degree m. Since q m (x) is a real polynomial, we 
write 
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q m {x) = a m x m + b m x m 1 + ... 

where a m # 0. For m = 0, (3) reduces to a 0 = cc no k 0 using the form for p n (x ) given in (1) 
and a 00 is uniquely determined. 

For the induction hypothesis, suppose that for in > 1, we can write any polynomial 
of degree m — 1 as a linear combination of { p n (x) : 



m— 1 

<7m-lW = Yj Hm-\ Pi( x )~ 

1=0 

Since q m (x) — ( ajk m )p m (x ) = is a polynomial of degree m — 1 the induction hy- 

pothesis implies there is a representation 



m— 1 

Pm( x ) 1 Pii x )' 

(=0 

Now set a m>m = ( ajk„) and the result (3) follows. [Ref. 4: p. 33] 

Using the theory of Fourier series developed earlier, we next determine the coeffi- 
cients a,. m explicitly. For / = 0, 1, ... , m, let c,. m = and let 

Then by construction, { <fi,(x) }™ 0 is an orthonormal set of polynomials. Writing (3) as 

m 

<7 mM = Y Cl ’ m M*)’ 

i = 0 

we see that the right-hand side can be interpreted as a (terminating) Fourier expansion 
of q m (x). Hence the results of Example 15 in Chapter II, Section B may be applied. In 
particular, by (7) in that section, the coefficients c, m are given by 

c i,m = < Pm . 4>i > = f q m (x) w( x ) dx. 

J a 

Changing back to the old variables, 
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a i,m = h i < <7 m - Pi > = h i f <7m(- v ) P/W W M dx - 

J n 



We are now in a position to show that the orthogonality property 



(4) 



f Pn( x )Pm( x ) w ( x ) dx = 0 ,m*n 

j n 



can be expressed equivalently as, 



( 5 ) 



'b 

p n {x) x m w{x)dx — 0 , m< n. 
v a 



To see (5), substitute the form (1) for p m (x ) into (4) where m<n. The linearity of the 
integral gives (5). On the other hand, since x m is a polynomial, we can write x m as a 
linear combination of the orthogonal polynomials, so (5) gives (4). Note that (5) implies 
each p„(x) is orthogonal to every polynomial of lower degree. [Ref. 4: pp. 33-34] 

B. THREE-TERM RECURRENCE RELATION 

The three-term recurrence relation is a useful result which holds for any three con- 
secutive orthogonal polynomials: 

(6) p n {x) = {A n x+ D n )p n _ x {x)~ C n p n _ 2 {x) ,n = 2,3,4,... 



where A n , B„ , and C n are constants given by 



A 



n 




>0 , D n = A 





C n = 



k n hfi— 

(k n -\) 



2 

2 



^-2 

K-\ 



>o. 



The recurrence relation is valid for n = 1 if p_ x = 0 with C, arbitrary. In this case, the 
formula for A„ also holds for rt= 1. (In the contrapositive form, this statement is a 
powerful tool for showing that a polynomial set is not orthogonal.) [Ref. 5: p. 234] 

To prove this, we begin by considering p n (x) — (kjk n . x ) x /Vi(- v )< a polynomial of de- 
gree no greater than (n — 1). We expand it in terms of the orthogonal polynomials 
v > a the technique of the previous section to obtain 
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( 7 ) 



n - 1 

Pni-v - - L . X p n _ x {x) = > a ;> _i aW- 



x n - 1 






The coefficients a i%n _ x arc determined by 



(S) =-h t - - JL - < xp n _ x {x ) , A<Jf) > , i < n 

K n-\ 



Because 



< x p n _ x (x) , p l (x) > = [ X p n _ x (x) Pi(x) w(x) dx = <p n _ x (x) ,x Pi(x)> , 

** a 

it follows that for i<>n — 3, .v p,{x) is of degree no greater than (n — 2). Hence 

< P n -\( x ) , x p^x) > = < x p n _ x (x) , Pi(x) > = 0 ,i<n- 3, 

since p n . x {x) is orthogonal to every polynomial of lesser degree. Thus the constants 
- . a n - 3 ^,-i are all zero, leaving <x n _ 2jt _ x and a*.,,,.,. With this knowledge, (7) 

becomes 



Pn( x ) - -r 1 - X Pn _ x {x) = a n _2,n- X Pn-l(x) + 

K n-\ 

Setting A n = kjk n _ , , D n = , and C„ = - then rearranging terms gives (6). 

To determine C„ explicitly, we write 

C„ = - = A n h n _ 2 <p n _ x (x ) ,p n _ 2 {x)> 



from (S). Since 



xp n _ 2 {x) = k n _ 2 x n 1 + • 



In- 2 
'<n - 1 

1 

V. 



[_ k n- l' Y " + ••' ] 



n—7 

Pn-\( x ) + YjPj.n-l Pj( x ) 
j = 0 



we can write 
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A n ~ 2 

C n = h n _ 2 — - 2 - < to , />„_ , to + ^ p Jn _ 2 pj{. x) > . 

”- 1 /=o 



Using the properties of the inner product and the orthogonality property of the 
polynomials, we obtain 



Qi 4—2 a 

*+ r 



to. />„-!(*)> “ 



4-2 



^-1 



44-2 4-2 
(/>•'„_ I ) 2 4-1 



as given in (6). [Ref. 6: p. 8] 

Into equation (6) we substitute the expanded forms of the polynomials 
p n (x ) . /Vito » P„- 2 to from (0 ^nd the constants A „ , C n from above. Equating coefTi- 
cients of jc'' - 1 gives B n . Since k„> 0 and h n > 0, it follows that A„ > 0 and C„ > 0 and the 
proof is complete. 

The nonnegativity of the constants A n and C„ is important for the converse of the 
result in (6). Favard showed that the existence of a three-term recurrence relation in the 
form of (6) implies that the polynomials of the set are orthogonal with respect to some 
weight function over some interval using Stieitjes integration [Ref. 7]. 

We observe that in order to generate the polynomials of an orthogonal class with 
the recurrence relation (6), we need the sequences of constants A n , B„ , and C„ together 
with two of three consecutive polynomials in the class. Other techniques provide what 
some authors call a pure recurrence relation requiring only two of three consecutive 
polynomials to define the class, because the constants as functions of n are contained 
explicitly in the recurrence relation. The recurrence relation derived in Chapter I, Sec- 
tion A.l for the Chebyshev polynomials is an example. 

C. CHRISTOFFEL-DARBOUX FORMULA 

The Christoffel-Darboux formula is an important identity which can be derived from 
the three-term recurrence relation. The identity is 



(9) 



Y^ h jpj( x ) P j(y) 



j = o 



4 k n Pn + 1 (- v K(j) - Pn( x )Pn + 1 ( >') 

4+. x ~y 



To prove (9), note that from (6) we have 

/’y+ito = (4+1* + b j+\)pM) ~ c j+\Pj-\i x ) 
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which we rearrange to give 



(10) xpj(x) = -j— p J+ , (x) - pj(x) + pj _ , (x) 

/l j + 1 ' v + 1 ' 7+1 

and similarly 

(U) JV/ J') = “J— / , y+i(>’) “ "T^ /’/>') + ~T^-Pj-\{y)- 

A j+ i w+i w+i 

These recurrence relations are valid for j = 0 if we set C, = 0 and p_ ,(x) = p_ { (y) = 0. 

Multiply (10) through by p^y) and (11) by p/.v), subtract the results, and then mul- 
tiply through by h ; to obtain 

h/kj _ 

hj (x - t) Pj(x) Pj{ y) = -7 — [ p j+ 1 (x)Pj( y) - p j+ , ( y)Pj(x ) J 

K j+\ 

+ J —j~ ~ [ Pj-\(x)Pj(y ) - Pj-i(y)Pj(x)]- 
Summing over j from 0 to n yields a telescoping series 
V 1 h n k n 

(x -y) )kjPj(x)Pj(y) = ~r [ p n +\(x)p n {y) - p n (x)p n+l (y ) J 

j=o *"+' 

from which the identity (9) follows. [Ref. 4: p. 39] 

Now subtract and add the quantity p n (x) p n + t {x) to the numerator of the right-hand 
side of (9) and let ^ tend to x to obtain a limiting case of the ChristofTel-Darboux for- 
mula: 



j ^ 

(12) Yj h j[ PM ) ] 2 = ~r JL [ Pn + 1 W Pn(x) - p'n(x) P n +\ M ]• 

j = 0 *"+> 

We will use (12) in the next section. 

D. ZEROS OF ORTHOGONAL POLYNOMIALS 

In Chapter I, we observed that the zeros of the Chebyshev polynomials {7" n (x)} are 
real, distinct, and lie in the interval ( — 1,1). This principle extends to any class of 
polynomials orthogonal on the real interval [a, 6]. 
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To see this, choose «>0 and suppose that p n (x) is of constant sign in [a. 6]. Then 
< p„(x ) , p 0 (.x) > ¥ 0, which contradicts the assumed orthogonality. Thus by continuity 
(and the Intermediate Value Theorem) there exists a zero x, e ( a,b ). 

Suppose that x, is a double root. Then p n (x)/(x — x,) 2 would be a polynomial of de- 
gree ( n — 2) and so 

o = < p n {x) , p n (x)l(x - x,) 2 > = < 1 , (p n (x)/(x - X!)) 2 > > 0 

which is a contradiction. Thus the zeros are simple. 

Now suppose that p n (x) has exactly j zeros x,, x 2 , ... , x, e ( a,b ). Then 

P n (x)(x - x,)(x - x 2 ) ... (x - Xj) = q n _j (x)(x - x,) 2 (x - x 2 ) 2 ... (x - x y ) 2 



where q n _j (x) does not change sign in (a, b) and 

< Pn( x ) . (x - x^lx - x 2 ) ... (x - Xj) > = < q n _j (x) , (x - x,) 2 (x - x 2 ) 2 ... (x - Xj f > . 

Since 

< Pn-J (- Y ) . (x ~ Xifix - x 2 ) 2 ... (x - Xj) 2 > ¥=0 

and 

< p n {x) , (x - x,)(x - x 2 ) ... (x - xj) > = 0 for j < n, 

then it must be that j > n. The Fundamental Theorem of Algebra precludes j > n and so 
we conclude j = n. [Ref. 5: p. 236] 

Thus all the zeros of p n (x) are real, simple, and lie in the interval ( a,b ), and so may 
be ordered 



a ^ ,n ^ x 2,n ^ ^ Xfc n < •" < X n n < b. 

The interlacing of zeros of /?„(x) and p„ + ,(x) follows from the Christoffel-Darboux for- 
mula. Recalling that the leading coefficient k n is positive for all />„(x), then (12) gives 

(13) p' +{ (x) p n (x) - p’ n (x) p n+l (x) > 0 , -oo < x < OO. 

Let u and v be adjacent zeros of p n {x). 'Then 

(14) p' n (u) p' n {v) < 0 
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since the zeros are simple. At these zeros, inequality (13) reduces to 



-/>'*(“) A ,+ i (»)>0 



and 



- /> ) P n+ ,(v) > 0. 



Multiply these two inequalities together with (14) to conclude 

/>„+,(«) /> n+1 (v)<0, 

so p n+l (x) has a zero between each pair of consecutive zeros of p n {x). 

Now let x„„ denote the largest zero of p n {x). Observing that p„(x) -* oo as .v -* oo, 
we must have p’ n (x„J > 0, and so by ( 13) 

Pn+\( X n,n) < °- 

But a +1 (jt) -» oo as x -* oo, so 1 (jc) must have a zero to the right of x nj> . Similarly, 

p n+l (x) must have a zero to the left of a: lin , the smallest zero of p n {x). Thus all n + 1 zeros 
of p n+l {x) are accounted for and interlace those of p n (x) [Ref. 8]: 

a ^ ^ x \,n ' Y 2^i+1 ^ "" ^ x k,n+\ x k,n • x /c+l,/i+l ^ ^ ^ x n,n -Tj+1,/7+1 ^ 

E. GENERATING FUNCTIONS 

The function F(x,t ) having a formal power series expansion in t 



is said to be a generating function for the set {/X*)} [Ref. 9: p. 1 29 j. For an appropriately 
chosen generating function, the generated set {/(.v)} is a class of orthogonal polynomials. 
By defining a class in this way, properties of the polynomials can be derived from the 
generating function itself. 

For example, consider 



OO 



r(x,t)-fuxy 



oo 



( 15 ) 
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which for fixed .v is the Taylor scries of F(x,t ) centered at r = 0. If we restrict .v to the 
interval [ — 1,1], then by considering the singular points of F(x,t ), we may conclude that 
the series is convergent for |t| < 1 [Ref. 4: p. 28]. Thus we can determine J'„(x) for 
n = 0, 1,2, ... by 

(16) / » w "i!rl7 r K | ~ 2 *' + ' ) I- 0 - 

Equation (16) yields f 0 (x) = 1 , f(x) = x, f 2 (x) = x 2 — , etc. We also note from 

(16) that 

- [” !( 1 
= 1 . 

We now derive some basic properties of the f n (x) from the generating function in (15). 

1. Recurrence Relation 

Differentiating (15), we find 

oo 

(17) -|£- = ,(l-2W)- 3 ' ! = £/,,■(*)<" 

n=\ 

oo 

(is) ^r ={x ~ 'K 1 - + 'T 3/2 = E 

n= 1 

Since (x — i ) — 1 ~r~ = 0, we have 

ox 01 



OO oo 

~ ~ 'E nfnW 1 " 1 = 0 

n = 1 n - 1 



which becomes 



E x fn'i x )t n - E N fn (x )‘ n = I J L'(xV +l . 

n = 1 n—\ ;i=l 
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Since f u {x) = I, we can start the sum on the right side at n = 0, then re-index so that it 
starts at n = 1 again. We then find 

oo oo 

'Yj [ x fn(x) - n U x ) ] t" = Yjfn-l M ‘ n - 

n=\ n—\ 

Equating coefficients of t" then gives 

(19) xf n '(x) - nf n (x) (x) , n > 1. 

Rewrite (17) as 



(20) ( 1 - 2xt + t 2 )- 3 ' 2 = * °- 

*=i 

Substituting the appropriate expressions from (20), (18), and (15), respectively, into the 
identity 

(l - f 2 )(l - 2 xt + t 2 )~ 312 - (2/) (x- r)(l - 2xt + t 2 )~ 312 = (l - 2xr + r 2 ) _I/2 . 
now gives 



(i w 

n - 1 



oo oo 

»=t n=0 



which, upon rearranging, becomes 



£/„+, w<" - E/„'-i - Z 

rt=l n = 1 /i=l n=l 



Equating coefficients of t" and gathering terms yields 

(21) (2n + 1) f n (x) =/„'+, (x) (x) , n > 1. 



Substituting (19) into (21) gives 

x fn’( x ) =/„'+, (-V) - (« + l)Z,(x) 
which by a shift of index front n-* n— 1 becomes 



30 



(22) 



*fn-\ (*) = fn'(x) ~ »/ n -|W . » > 2. 



Substituting again Prom (19) gives 

(23) (x 2 - l)f n '(x) = nxf n (x ) - nf n _ ,(x). 

Multiplying (19) through by (x 2 — 1), we have 

.t(.v 2 - 1 )f n '(x) - n(x 2 - 1 )f n (x) = (x 2 - 1 (x). 

Now substitute for (x 2 — 1 )/,'(*) and (x 2 — 1 )/,'_, (x) from (23) and gather terms to get 

(24) nf n (x) = (2 n - 1) xf n _ x {x) - (n - l)/„_ 2 (x) , n> 2. 

Equation (24) is a three-term recurrence relation for {/X-t)}. The advantage of 
this form is that beginning with f a (x) and y[(x), we can now generate any member of 
(X(x)} by iterating (24) and thus avoid the differentiation in (16). Note that since 
f 0 (x) = 1 , yj(.r) = x, (24) implies that {/^(x)} is a set of polynomials. [Ref. 9: pp. 159-160] 

2. Ordinary Differential Equation 

We continue the same line of reasoning to extract additional information about 
U(x)}. Differentiating (22) yields 

(25) x/;_, (x) =f'\x) - (n + 1 )/„'_, (x). 

From (19) we have f ‘_ , (x) and, after differentiating, (x). Substituting these ex- 
pressions into (25) gives 

x [xfn’{x) -(n- l)/„'(x)] =f n "(x) ~(n+ l)[x/„'(x) - «/,(x)] 
which when rearranged becomes a second order ordinary differential equation 

(26) (1 - x 2 )f n "(x) - 2xf n '(x) + n(n + 1) f n (x) = 0. 

The {/ n (x)| are solutions of (26) for // = 0, 1,2,.... [Ref. 9: pp. 160-161] 

3. Orthogonality 
Rewriting (26) as 

(27) [(1 - x 2 )/.'(x)] + n(n + 1 )/„(x) = 0 
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we recognize the structure of a Sturm- Liouville eigenvalue problem. Since the points 
.r = ± 1 are singular points, we require that/fcc) andy^'(.v) be finite as .t -*• ± 1. From the 
associated theory of the singular Sturm- Liouville problem, we conclude that the \f„(x)} 
are orthogonal on the interval [ — 1,1] with weight function vv(x) = 1. 

Alternately, we combine (26) with 

[0 — -* 2 )/m'W] + m ( m + = 0 



where n ^ m, to obtain 

(25) frn(x) [(1 - -^ 2 )/n'W] ~ to) [0 ~ * 2 )/m'W] 

+ [n(n + 1) - m{m + l)]/„(^)/ m (.t) = 0. 

Since 

-£■ [0 ~ ^K/mW/n'W -/m'W/nW}] 

= -^- [/mW{( 1 -* 2 )fn'(*)} -fn ( x ){( 1 -^ 2 )/m'W}] 

-/»(*) [(1 - * 2 )/«'to] ~to) ■£-[(!- -v 2 )/ m '(-v)], 



we can write (28) as 

■£•[(1 - -U W/,W)] + [n 2 - nf + n - rn]/„(.t)/ m W - 0 



or 

(n - m){n + m+ 1 )f n {x)f m {x) = -^ [(l - x 2 )[f m '{x)f n {x) -/„,(*)/,/(*)}]. 
Integrating from x = — 1 to x = 1, we obtain 

(n - m){n + m+ 1)J f n (x)fjx)dx = [(l - Jf 2 ){/„'(*)./!»(*) 

Since (1 — .t 2 ) = 0 at jc = ± 1, we have 
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(h — m)(n + m + 




Recalling that n =£ m, we conclude 



f /nW/mW* = 0 

J -i 

i.e., the polynomials {^(.v)} are orthogonal on the interval [—1,1] with respect to the 
weight function w(x) = 1. [Ref. 9: pp. 173-174] 

These are the Legendre polynomials that were previously defined using the 
Gram-Schmidt process. Thus the f„(x) of (15) are in fact P n { x) and (15) may be written 



(29) (1-2 xt + tY l7 = Y J P nix)t n , 

n= 0 



establishing the equivalence of the generating function definition with the Gram-Schmidt 
definition. 

Legendre and Laplace concluded that the P n (x) in (29) were polynomials of de- 
gree n in the variable x by examining a series expansion of the function 



[«/ 2 ] 






n=0 



m = 0 



They reasoned the orthogonality directly. From (29) we write 



(l— 2 xr+r 2 ^ ^ Pj(x) r 

(=0 



and 



(l— 2xs+s 2 ) '^ 2 = ^ Pj(x) .</. 

7=0 

Multiplying these power series together via the Cauchy product formula yields 
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1— 2xr+r 2 N /l— 2 xs+s 2 n=o k=o 



on ri 

-EE w <■* *"■* 



-EE w w 



;i=0 L /+y=n 



rV. 



Integrating from -1 to 1 with respect to x gives 



f -= 






2xr + r 2 >y 1 — 2x5 + 5 2 «= o /+/=« 



oo r p i 

-EE J ww 

n-i) Mm ~ 1 



dx 



rV. 



Through tedious calculation, the left-hand side of this expression becomes 



-Jrs 



i , i + 

log' 



1 - Jrs 



a function of the product rs. From this we conclude that the coefficients of the terms 
in the series on the right-hand side are zero when i and j differ, i.e., 



r pix) Pj(x) dx =o , 

J -1 



and the orthogonality is established. Finally, 



(i-w)-^|„,.- i J T -E'" 

n = 0 

or /*„(!) = 1 and so the P n (x) are in fact Legendre polynomials. [Ref. 8] 

In a similar fashion, the norm of the Legendre polynomials can also be obtained 
from the generating function. Let 



= K ‘ = f [ W f dx. 

J -i 

Square both sides of (29) and integrate with respect to x on [ — 1,1] : 
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! 




P m (x) P„(x)dx 




1 — 2xt + t 



dx. 



Dy orthogonality, the left-hand side is zero unless m = n, while the right-hand side is 
integrablc in closed form: 



In (l — 2xi + < 2 ) | ^ 

This function can be expressed as a difference of two logarithms, each of which has a 
convergent Maclaurin scries expansion in ( —1,1). When combined this yields 






/?=0 



z 



n = 0 



C n t 



2n 




2 

2 n+ 1 



) 



t 



In 



Comparing coefficients of t 2 " on both sides gives c„ = — . 

Since there is no systematic theory for determining generating functions, finding 
one for a polynomial class can be a problem. The work above bears this out. Unfor- 
tunately, the proofs above do not easily generalize to related classes. With this in mind, 
let us summarize the key steps in proving that (15) is a generating function for the 
Legendre polynomials. First we established that f„{ 1) = 1 for n > 0. Next, we showed 
that 



(l ~ 2ar/ + f 2 ) 1/2 = YjfnU) i n 

n = o 

for f n {x) a polynomial of degree n in the variable x. Finally, we showed that these 
polynomials were orthogonal on the interval [ — 1,1] 

f f,(x)f m ix)dx = 0, m t- n. 

J -\ 

These three points are sufficient to show that the generated class is the Legendre class 
of orthogonal polynomials. 
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We now provide an alternate proof attributed to 1 lermitc. Our interest in this 
proof is mainly the technique which suggests a method of generalization that we will 
take advantage of in the next chapter. 

We begin with ( 15) reproduced below 



(l-2w + > J )- ,p 

n=0 

Multiply both sides of this equation by x*and integrate from -1 to 1 to obtain 



(30) 



f 

J -1 



x k dx 



77^ 



XI + l 



2 



x k / n (x) dx. 



n = 0 _1 



Now change variables from jc to j/ via 

(l - 2 xt + t 2 )~ U2 = 1 - ty 



giving 



JC = 



'(i-V) 



+ y, dx = (\ — ty)dy. 



The left-hand side of (30) becomes 

[t '(■-/) 



I 

J -i 



+ T 



(1 - v) 



(1 - ty)dy 



or simply 



r, [i '(' -/)+>•]* 4 --. 



Expanding the integrand of this last expression via the binomial Theorem, we obtain 



k 



,i (* ->■)' k-j . 

jj'J *■ 



y=0 
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When written as 



k 




J 



we can identify the form as a polynomial of degree k in the variable t. Comparing this 
form to the right-hand side of (30), we conclude that 



i. Q.,f n (x ) is orthogonal to all polynomials of lower degree. The same argument as before 
gives^(l) =1 for n > 0 and the proof is complete. 

F. HYPERGEOMETRIC SERIES 

The term "hypergeometric" was used in 1655 to distinguish a series that was "be- 
yond" the ordinary geometric series 1 + x + x 2 + ••• . In 1812, Gauss presented the power 
scries 

ab x 1)W>+ 1) ,v 2 a(a+l)(a + 2)b(b+l)(b + 2) , v 3 

+ c 1 ! + c(c+l) 2! + c(c+l)(c + 2) 3! + 

c # 0,-1, —2, ... which is known as Gauss' series or the ordinary hypergeomeiric series 
[Ref. 10]. 

Convergence of this series for |.v| < 1 follows directly from the Ratio test. By 
Raabe's test, convergence can be shown for |.r| = 1 when (c — a — b) > 0 [Ref. 1 1: p. 5J. 
Gauss also introduced the notation 2 F x [_a,b ; c ; .t] for this series. Note that 
2 F,[a,6 ; c ; jc] may be considered as much a function of four variables as a series in .r. 
[Ref. 12: p. 1] 

With the shifted factorial, the ordinary hypergeometric scries can be expressed 



x k f n (x) dx = 0 , for n> k. 



-l 



oo 
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Below are some examples of important functions which can be expressed as ordinary 
hypergeometric series. 

Example 1: log(l + x) = x 2 F',[l, 1 ; 2 ; — x] 

Example 2: sirr'(x) = x 2 F,[ 1/2, 1/2 ; 3/2 ; x 2 ] 

Example 3: tair'(x) = x 2 F,[l/2, 1 ; 3/2 ; — x 2 ] 

The generalized hypergeometric series is formed by extending the number of pa- 
rameters, an idea attributed to Clausen [Ref. 1 1: p. 40]. 



r F$- [T l ? ^2’ > ^1> ^2 



oo 

, ... , b s ; x] = 



(^|)ai(^2)/i ••• ( a r)n 

(b 1 ) n (h 2 ) n ...(b s ) n n\ 



n = 0 



Note that since (a) n ^J(a) n = n + a, a hypergeometric series Z c„ x" is characterized by the 

n-0 

fact that the ratio c nJ .Jc n of coefficients is a quotient of two polynomials in the index n, 
i.e., a rational function of n. 

The Ratio test can be used to show convergence for all values of x when r <; s and 
for | x | < 1 when r = s+ 1. When r > s + 1, the series diverges for all x^O and the 
function is defined only if the series terminates. The series terminates when one or more 
of the numerator parameters u, is zero or a negative integer [Ref. 1 1: p. 45]. This is an 
important characteristic of the hypergeometric series that will be used later. A power 
series that terminates gives a polynomial which is defined for all x. In this case, the pa- 
rameters />,, ... , b , may be negative integers as long as the series terminates before a zero 
is introduced into a denominator term. 

Examples of the generalized hypergeometric series include familiar functions such 
as: 

Example 4: (1 + x) a = ,F 0 [ —a ; — ; — x] 

Example 5: e* = 0 F 0 [ — ; — ; x] 

Example 6: sin x = x 0 F,[ — ; 3/2 ; — x 2 /4] 

Example 7: cos x = 0 F,[ — ; 1 /2 ; — x 2 /4] 

Example 8: The Bessel function of order a 

W2)V,[- ;«+ I; -a: 2 /4] 

J - (x) nTTn “ 



where dashes indicate the absence of parameters, i.e., when r = 0 or s = 0. We will also 
use a common alternate notation 
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A 



A*- 



>«r 




X 



for either the ordinary or generalized hvpcrgeometric scries. [Ref. 12: p. 4] 
1. Chu- Vandermonde Suin 
The Chu-Vandermonde sum 



2 f,[ -n, a ; c ; 1] = 



(c ~ a )n 

(c) n 



is one of many useful summation formulas. Since this one will be used in a later chapter, 
the proof is provided below. 

Basically, this is a consequence of the General Binomial theorem 



OO 




k=Q 



Starting with the identity (1 — .v) _ °(l — x)'* = (1 — expand both sides. Us- 
ing the Cauchy product on the left side and the General Binomial theorem on the right, 
we have 



oo oo . 

V n V ( a + b)„ „ 
c„x - ? — .V , 



n=0 



ho ^ 



where c„ = £ ^ - . Equating coefficients of x n , we have 

*=0 (!)*(' )n-k 



y WMn-k _ (a + b) n 
{y (•)/:( On-Zc (On 



In order to express the left side as an ordinary hypergcomctric series, 2 F U mul- 
tiply both sides by (1)„ and use the identity (1)J( 0-.-* = ( — 0*( — «)* to obtain 



( 31 ) 



I 

A'=n 



( -n) k {a) k 
( 0 * 



(-0 fc W„-* = (« + *)„ 



Next, 
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(-l)Wji -(-!)"( 

= I) ...(* + »-*- I) 

= ( -1)"( -b)( -b - 1) ... ( -b - n + k + 1) 

_ ( ,y. On 

1 } (-b-n+l) k • 

Using the above result, equation (31) becomes 



I 



(-")k( a )k 1 n « (a + b)„ 

(1 )/< (-b-n+ l) k 1 j (-b-n+ 1)„ 



Let c = — b — n + 1, and substitute to get 



2 F\l~n, a;c; 1] = (-!)" 



( — c + a — n + 1),, 



(—!)"( -c + a — n - f 1) ... ( — c + a) 
ton 



(c — a + n — 1) ... (c — a) 
ton 

(c - a)n 
ton ' 



n + Ufe 

n + 1)^ 
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IV. JACOBI POLYNOMIALS AND SPECIAL CASES 



This chapter focuses on the orthogonal class known as Jacobi polynomials, P^' p \x). 
Until recently, the classes of orthogonal polynomials considered "classical" were usually 
given to be Jacobi, Gcgenbaucr (also called ultrasphcrical), Chebyshcv (of first and sec- 
ond kind), Legendre (also called spherical), Laguerre, and Hermitc. The Jacobi 
polynomials hold a key position in this list since the remaining classes can be viewed as 
special or limiting cases of this class. Today, the classical orthogonal polynomials arc 
taken to be special or limiting cases of cither of two very general orthogonal classes 
known as the Askey-lVilson polynomials and the q- Racak polynomials, between which 
we will establish a formal equivalence. Because of their complexity, these classes arc 
described in Chapter VI after the necessary additional theory has been developed. 

The results derived in the text that follows arc arranged in tabular form by class at 
the end of this chapter. 

A. JACOBI POLYNOMIALS 
1. Definition / Orthogonality 

The Jacobi polynomials P^' p \x) arc generated by applying the orthogonalization 
step of the Gram-Schmidt process to the standard basis (1, jc, x 2 , ... } of P[ — 1,1], with 
respect to the weight function given by a continuous beta distribution on [ — 1,1] 

w(x;a, (1) = (\ - -rf ( 1 + xf 



for a > — 1, /? > — 1, i.e., 



f 1 p { :'%)(i-xn\+xfjx=[h ( ;- /J) ]- 1 & mA . 

j ~\ 



The Jacobi polynomials can also be represented by hypcrgcomctric scries 



(i) 



it%)- ( “ + l) " 



n\ 



i 2 ‘ i 



E, 



-n, n + <x + /? + 1 1 — x 

a + 1 '~T~ 



This set of polynomials is thus standardized (as was done for the Legendre class in 
Chapter II, Section B) : 
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We shall demonstrate that these two characterizations of the Jacobi polynomials arc 
indeed equivalent. 

To show that the polynomials defined in (1) arc in fact orthogonal with respect 
to vv(.v ; a, /?) = (1 — Jt)’(l + x) p on [ — 1,1], it suffices to show that P^' p \x) is orthogonal 
to one polynomial of degree m for m = 0, 1, 2, ... , n — 1. This is because any P^' p \x), 
0 <1 m < n — 1, can be expressed as a linear combination of such polynomials. While any 
polynomial of degree m could be used, we choose (1 + *) m for reasons that will become 
apparent. 

To establish orthogonality, we consider 



< P„’ ^\x ) , ( 1 + .x) m > = Pf’ p> (x) ( 1 + x) m ( 1 — x)* ( 1 + x) p dx. 

“'-I 



:>(*. P)( 



By (1), the right-hand side becomes 



(«+0n V 

n\ l_j 



k=0 



( ~n) k (fl + a + /? + l) fe 
(«■+!)»*! 2* 



f'd-.rf+’d 



+ x) m+ ^ dx 



Using the last result of Chapter II, Section C.l, the change of variable „r = 1 — 2t yields 



(« + 1), V 

n\ l_j 



k = 0 



( ~^)fe ('? + « + /? + 1)^ \ 

(a+l) k k'.2 k ) 



D(k + a + 1 , m + (t + 1) 



( a + 1 )/i 



n\ 



f t 

E 

k = o 



( — (ft + a + /? + l) fc \ ^ m+0L+ p + x V(k + a + 1) r(m + /?+!) 

k\ J (a + l) fe f(/: + m + a. + (] + 2) 



Identities for the gamma function allow us to simplify this expression to the form 
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2*+(J+m+i r(a + n + i) rfm + (t+l) 

r(m + a + p + 2) 



l 



( —n) k (n + a + /? + Oft 
(m + a + ft + 2) k k\ 



which can then be written 

—n , n + a + /? 4- 1 

; i 

m + a + P + 2 



2 * +p+m+1 r(a + /i + 1 ) f(m + /?+!) 

r(m + a + /? + 2) 



The Chu-Yandermonde sum from Chapter II allows us to write this as 

2 *+p+m+i r(g + n + [) f ( w + /?+!) (m + 1 - n) n 

r(m + cc + /? + 2) (m + ce + fl + 2) n 



Thus, 



< P^' p \x),{ \+x) m > 



jc n , m = n 

0 , m = 0, l, 2, ... , n — 1 



where 



n\ 2 n+ * +p+l r(« + a + 1) Un + p + 1) 

° n ~ Y(2n + a + p + 2) 

which justifies the orthogonality. It is possible to extract the value of by modi- 

fying this argument, but we defer this computation till the next section, where it will be 
easier. 

2. Ordinary Differential Equation / Rodrigues' Formula / Norm 

We begin deriving the ordinary differential equation for the Jacobi polynomials 
by noting the general formal result: 



( 2 ) 



_d_ 

dx 



, . . 


• , a r 


a, . 


•• a r 


a \ + 1 , . 


» a r + 1 




i ' X 

•A J 


-v 


b r S 

.. u s 


J>\ + 1 1 •• 


; x 

• , h 5 + 1 



This can be seen from the definition of the generalized hypergeometric series by 

writing 



43 



_d_ 

dx 



fs 



«l, ••• . a T 

J> !•- - b S 



_d_ y («i) fc - (g r ) fe ,y fe 

^ ^ (^i)a — (^r)/E k' 

y Me - fa)* M 

tt Mk-(bs) k (k-l )! 

y fa)*+l -fa)*+l ** 
^ fa)/<+l — fa)*+l k\ 



Noting that (a)* +l = a (a + 1)*, we have 



_d_ 

dx 



A 



a \y ••• > 

J>\, ••• , b s 



X 



<*1 y fa + l)*-fa + 1) A y 

b\--b s ^ fa + 0* ••• fa + 1) A k\ 



from which (2) follows. Differentiation of the hypergcomctric series is justified by re- 
calling that the power series is convergent for [ jc| <1 when r = s + 1. 

We apply this result to the Jacobi polynomials to obtain 



_d_ 

dx 



( (a + 1]* 



P^'P)( X \ = -2— 
n W dx 



1 L 1 






— n, n + a + /? + 1 
a + 1 

(a + 1)„ — n (n + a + /? + 1) 



n\ a + 1 

(n + a + P + 1) (a + 2) n _ { ( n — 1)! 



2 (n — 1)! 



(a + 2)„_, 




2 r 1 



—n + l,tf + oc + /? + 2 
(x + 2 



p^’ p+ %) 



1 — .V 
2 



which simplifies to 



( 3 ) 



jL n(«./?) r .A (n + ct + P+l) 
dx n [ ) ~ 2 



^iV’ /?+1) w 



With these results established, we now consider the orthogonality property of 
the class: set 



J -1 

By (3), we can write 
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*'-1 

Integration by parts yields 

< 5 > w-zfhr /.',[*■ (^ n M(i — )■(• 

noting that the boundary terms vanish, so that 

( 6 ) 4^ = ' w+ ~ 2 + J J ( p m+l’ P ~'\ x ) %+\( x ) 0 ~ x) a ~ l (1 + x'f~ ] dx 

where 

( 7 ) %+ lW = 0 -^ 2 ) P ( n' P \x) - a (1 + *) P ( n' P \x) + 0(1-*) Pn' P \x) 

is a polynomial of degree ( n + 1). Thus by (6) we may express < 7 „ + ,(jc) as a linear combi- 
nation of Jacobi polynomials 



n + 1 

%+\ W = Yj c j P j t ~ x ' 

;=o 

We would now like to show that c 0 = c, = ••• = c„ = 0, so that only the last term 
of the summation survives. The constants c ; , j = 0, 1, 2, ... ,n can be determined by 
substituting this expression for q n ^{x) into (6) and using (4). 



/(*. P) 

1 m, n 



■' /><*. <»(.<•) />£' P\x) ( I - jr)* ( 1 + Jr)'* dx 



-2 



^ r* 



m + a + fl £-, Cj ) , ' m+1 
y=o ~ l 



J' 



(a-l, pta-l, /?-!) 



For each m<n— 1, by (4) the left-hand side / m , 
zero unless y = m -f 1 <; n. Thus 



(*)(1 -.^“'(l 4-^‘Vv. 

, = 0, but the right-hand side is 
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which means from (6) 



w c m + 1 L'w+l J ’ 



cj = 0 , j< n 



( 8 ) 



/(*./?)_ 1 1 ) 
m < n *+1 ‘m+X.n+X 



and 

( 9 ) 



<?„+ iW = C «+1 ^+/' p ]) W- 



To determine the remaining constant c n+ „ we note that by (7) and (9) 
(l-.r 2 ) -£■ P^ p \x) - a (1+x) /f >/7) (x) + /?(!-*) /f >/?) (x) = c„ +1 PfJ' P ~%). 



Letting x= 1, then from the hypergeometric series definition of the Jacobi polynomials 



we have /**’ w (l) = — — which gives 



n\ 



-2a 



( a +0* 



I C n+ X 



(“)« 



+ 1 



«! 



(«+!)! 



and so 

( 10 ) 



c n+ x = -2(n+l). 



Combining (5), (6), (9), and (10), we have 

(11) -£■ [ Pt P \x) (l-xr (l+xf ] = -2 («+l) [ pfc?' P ~%) (1-x) 3 - 1 (1+x)^" 1 ]. 

To obtain the second order ordinary differential equation, change n -* n— 1, a->a+l, 
and /?-*•/?+ 1 in ( 1 1 ) to give 

J_ |>(*_+i. P+ D W (l_x)* +1 (1+x)^ 1 ] = -2 n Pt' P) (- x) (l-xf(l+xf. 



Then by (3), 
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d 2 _d_ 

d.X ll+OL+P+l dx 

i. e.,y = Pn' p \- v) satisfies 



/f^.v)(l— .vf +1 (l+.v/ +1 



= -2 >i P ( n' P \x) ( 1 -xf ( 1 + x f , 



_d_ 

dx 



(l-xf +l (1+4 



,p+l 



dy_ 

dx 



—n (n+cc+P+ 1) (1— 4* (1+4^+- 



Using the product rule to expand the left-hand side, we have 

d-44i + 4^ («+i) (1+4/ + (/?+ 1) (1-4/ + (i--* 2 )/' ] 

= — n (n+a+p+l) (1—4’ (1+4^ + 

which, when simplified, becomes the second order ordinary differential equation for the 
Jacobi polynomials y = P^' p \ x) : 

(12) (l-* 2 )/ + [ (/?“«) - («+/?+ 2) x]y’ + n ( n+a+p+l)y = 0. 

The reader is invited to compare this result with the second order ordinary differential 
equation for the Chebyshev polynomials {T n (4) given by Equation (5) in Chapter I. 

By iterating ( 1 1) k times, we obtain 

"4- [ P ( n' P \x) (1-4’ (1+xf ] = (-D* 2 k (# 1 + 1 )* [ />£?• P ~ k \x) {\-xf- k (1 +xf~ k ]. 
dx 

Setting a -* a +k, and /? -*• p+k gives 



(l-4“ (1+4* Pn +k p \x) = (J T [ (l-4’ + * (i+x) p+k Pn +k - p+k \x) ], 

7 k {n+\) k dx K L J 

or equivalently in terms of the weight function, 

4-V ; a, P) Pn+ P \x) = !~ U< -4 [ VV(X ’ a + k ' £+*) p n +k ’ P+k \x) ]. 

2 (/i+l)^ dx 



the general Rodrigues' formula for the Jacobi polynomials. Letting n = 0 and noting 
Pf p \x) = 1 gives 



( i —4“ ( i +x) p p k ’ p \x) = 4^~ 4r [d-4 1+fe (i+4 /?+fe ]. 



k\ dx' 



47 



or 



vv(.x ; a, P) P ( f %) = ^-jp- [ w(.\r ; u+k, p+k) ], 

2 k\ dx 

the classical Rodrigues' formula for the Jacobi polynomials. 

Finally, we can use these ideas to obtain the value of By direct com- 

putation, 



[/,<*’ 0>] 1 = f (1 -xf ( 1 +xf dx 
J -i 

2 a+/?+1 r(a+i)r(/?+i) 

V(a+p+2) 



We will use this below. Combining (8), (10), and (4) when m = n yields 




n+a+p L Wl J 



Making the changes n — *• n— 1, a -*■ a+1, and /? -> /?+ 1, this may be rewritten as 






Iterating this relation n times produces 



W-T 



(n+a+P+ 1)„ r , (a+n, 0+n)~]—l 

„i n . L /! 0 J 

4 n\ 

r( 2 »+a+/?+l) 2 2n+ * +/?+l r(//+a+l) r(/;+/?+l) 

2 ln n\Y{n+a+p+\) V{2n+a+p+2) 

2* +/?+1 r(/t+q+l)r(/;+/?+l) 

(2«+a+/?+l) n! r{n+a.+P+\) 



3. Generating Function 

As was mentioned in Chapter III, finding a generating function for a class of 
orthogonal polynomials can be a challenging task. Fortunately, a generating function 
for the Jacobi polynomials can be found by mimicking a technique attributed to Ilermite 
in his work with the Legendre polynomials. (See Chapter III. Section H.3.) 

We seek a generating function 
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(13) 



F(r, x) = 'Y J C n P ( *' ti Xx) r n 



n = 0 



for constants c„ , n — 0, 1, 2, ... , with 



f Pjj 3 ’ ^(x) x fc ( 1 — x)* ( I + x'f dx = 0 , k < n. 

3-i 



To begin, we multiply (13) by x k (1 — x)’ (1 + x) p and integrate from -1 to 1 to obtain 



M ^ r M 

(14)1 x k F (r,x) ( I — xf ( 1 + xf dx = V' c n f P^’ ^(x) x k ( 1 — x)“ ( 1 + x'f dx 
J -i >i=o ^-i 



Note that the summation of the right-hand side is from 0 to k because the orthogonality 
property makes each term zero for n> k. We observe that the unknowns to be deter- 
mined in (14) arc F{r,x) and c„. When one is given, the other can be found, so we first 
consider the left-hand side of (14) 



f x k F (r,x) ( 1 — x)* ( 1 + xf dx. 
3-i 



Setting J 1—2 xr+r 2 = l—ry and substituting for x = y + ^ — —— j r yields 



I 



■ r , 0 V) 

-,L y+ — 



F[ r, y+ 



(IV) 



i -y- 



OV) 



\+y+ 



(1 -/) 



2\ 



( 1 -ry) dy. 



Factoring out the terms (l— y)’ and (1+y)^, we have 



(15) 



-i L 



v+ 



0 -/) 



2\ -n 



(l-yr(l +yf F I r,y+ 



(1 V) 

2 

( I +>■) r 



1 + 



( 1 -y) r 



( 1 -ry) dy. 



We would now like to make a judicious choice of F in order to facilitate calcu- 
lations. Following the lead of the Legendre polynomials (Chapter III, Section L.3) 
suppose F is such that 
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(16) 



l : [ r,y+ 



o-y 2 ) 



(1+y) r 



1 + 



(1 -y)r 



(l-o0 = 1, 



then the integral in (15) becomes 



f 

J -i 



T+ 



(>-/) 



(1 -y)*{\+yf dy. 



We note that the integrand is a polynomial of degree k in the variable r as is 
the right-hand side of (14). From (16), vve get 



/ ( 1 ->’ 2 ) r \ 

F[r,y+ K ■/ = (1 -ry) 



( 1 r 



1 + 



0 -y)r 



-p 



which becomes, via l— 2xr+r 2 = (1— ry) 2 , the generating function for the Jacobi 
polynomials: 

F(r,x ) = 2 X+ ^ (l— 2 .vr+r 2 ) -1 ^ 2 [l— r+(l— 2.rr+r 2 )^ 2 ] - * [l+r+(l — 2.xr+/' 2 ) 1 ^ 2 ] - ^. 
Using this generating function in (13) gives 



c n p£-%) r n = 2 a+ ^ (l— 2 jcr+r 2 ) ! ^ 2 [1— r+(l — 2.Y/'-t-r 2 ) 1 ^ 2 ] ^ [l4-r4-(l — 2-vr-Fr 2 ) lyf2 ] ^ 

( a +l) n 



*=o 



from which we can determine c„ by setting x = 1 and recalling that w ( 1 ) = 
Thus, 



m 



oo 



n i , \— 1 r->— «/• \— on P it \— a— 1 

r = 2 0-') l 2 (I-O J2 F = (l-r) 



> i=0 



which by the Binomial Theorem becomes 



(«+l )n n V ( a+1 )n „ 

r = / ; r 



V b*+U„ „ V 

L c *— t 7 ~ r = L 



i 1=0 



>?=0 



//! 



and so c„ = 1 for all /?. [Ref. 8] 
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D. SPECIAL AND LIMITING CASES 



With the structure of the Jacobi polynomials established, we turn now to the role 
of the parameters a and /?. 

1. Special Cases 

For certain choices of the parameters a and /?, we find that the classes previously 
examined and several new classes are produced as special cases of the Jacobi class. 
These subclasses inherit the structure of the parent class which often provides a direct 
way to establish specific properties (i.e., polynomial nature, orthogonality, etc.) 

Table 1 provides the choice of parameters for selected classes. 



Table 1. PARAMETERS FOR SELECTED CLASSES 



Class 


Parameters 


Jacobi 


a > — 1, /?>— 1 


Gegenbauer 


& 

II 

II 

V. 

73 


Chebyshev, 
First Kind 


« = /? = — 1/2 


Chebyshev, 
Second Kind 


II 

II 


Legendre 


O 

II 

II 

8 



2. Limiting Cases 

In this section, we briefly examine the Laguerre polynomials and the Mcrmite 
polynomials. Using the hypergeometric series definition of the Jacobi polynomials, we 
show how the Laguerre class is a limiting case of the Jacobi class. 
a. Laguerre Polynomials 

The Laguerre polynomials, L^\x), are defined in terms of hypergeometric 

series 



( 17 ) 




fn+a\ 


—n 




4 \x 


\ n r 1 


a+1 



The first relationship to establish is 
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lira P { n' P) 

p - OO 




To show this, we begin with the hypergeometric series representation for the Jacobi 
polynomials 



nf*. P) 

£ n 



{y) = 



n+CL \ 

24 

n ) 



— n, /i+a+/?+ 1 1 — y 

a+1 : ~ 



2 x 

When y = 1 — ~ , we have 



p(*. P) 

r n 





— n, rt+a+/?+l 
ft-) - 1 



x_ 

P 



Writing out the power series, we obtain 




In the limit as P -* oo, 




ratio 



(^+ a +A+l)/c x k 

(a+l)**! p k 

Q;+a+/?+ l) /c 
(a+1)* A! p k 



(n +a+fl+l) k 

P k 



and so 



lim Pj?' p) 

p-OO 




tgh 

(a+l) fc k ! 



x k -L { *{x) 



completing the proof. [Ref. 13: p. 103 [ 

To obtain the orthogonality relation for the Laguerre polynomials, we start 
with the orthogonality relation for the Jacobi polynomials 



I'm p \ y) p n' p) ( y) ( 1 -yY ( 1 +yf 4 = o , m*n. 
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Letting j/ = 1 — gives 



f 



?<*. p) 



P 



■L \ n(». P) 



Lx 



2.x \* 

P 



O-f) 



2 . ft 
— a„v = 0 , m ^ 



Passing the constants outside of the integral and dividing through by them leaves us with 



r ^ "O' - ■ ^ ) * '’( i_ ■ ) ** (' - 1 y * ■ - ° ■ " # *• 

Taking the limit as p -* oo, we obtain the orthogonality relation for the Laguerre 
polynomials 



(18) 



f 



L^\.x) L^\x) x* e x dx = 0 , m ^ n. 



b. Hermite Polynomials 

The Hermite polynomials, H„{x), are defined 



(19) 



//„(*) = (2*) Vo 



-n/2 , (-«+l)/2 i 



x 



In a fashion similar to that for the Laguerre polynomials, the Hermite polynomials are 
a limiting case of the Jacobi polynomials via the Gegenbauer polynomials. Specifically, 

II n {x) = n\ lim X~ n ^ 2 C^(}.~ l l 2 . r) 

A—*oo 

which allows a derivation of the orthogonality relation 



f°° r 2 

(20) H m (x) II n (x) e~ x dx = 0 , m * n 

^-oo 

from that of the Jacobi class. (Ref. 13: p. 107] 

C. DISCRETE EXTENSIONS 

We turn now to orthogonal polynomial classes which use the discrete inner product 
introduced in Chapter II, Section A. 2 
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<v 

< /ig > = y^/(-v) g(jc) w(.v). 

jc=0 



Here instead of a continuum the support of the weight function is concentrated of a fi- 
nite set of discrete mass points {0, 1, ... , N). The polynomial classes of particular in- 
terest are the Hahn, dual Hahn, and the Racah polynomials. 

1. Hahn Polynomials 

The Hahn polynomials - actually discovered by Chebyshev - were independently 
realized by physicists working in angular momentum theory via 3-j, or Clcbsch-Gordon 
coefficients [Ref. 8]. We define this class by the generalized hypcrgeometric scries 



( 21 ) 



Q n ( x ; P, 



n) = 3 f 2 



—n, —x, n+a+P+l 
a+1, —N 



1 



for a > — 1 , P > — 1, where N is a positive integer and n = 0, 1, ... , N. From the power 
series 



nr a ^ V (-n)k(-x) k (>t+*+P+ 0 * 

Q„(x , a. /J, 2, («+I) ,(— ,V),« () ’ 



we note that the variable x does not appear where we have come to expect. Since 



(-■*)*"(-*) (-■*+!) - (-v+A-l) 

= (-l) k (x)(x-l) ... (x-k+l), 

we conclude that Q n (x ; a, p, N) is a polynomial of degree n in the variable x. Because 
0 < n < N, this set of iV+1 orthogonal polynomials is finite for fixed a and /?. [Ref. 14: 
p. 33] 

The Hahn polynomials satisfy the discrete orthogonality relation 

v 

(22) ^jQ m {x\o., P, iV)Q„(x; a, /?, iV)w(x;oc, P, N) = 0, nv£n 

.x=0 

where the weight function is given by a hypergeometric distribution 
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w(jr;a, /?, A r ) = 



( 23 ) 



x+a\ ( iX—x+f 



(a+l) x (/?+l) v _ x 



(1), 



(D.v-, 



We note that the weight function may also be written 



(24) 



w(jc ; a, /?, /V) = 



(A+l),v ( a +0x )x 



U),v (I) x (-A-/?), 



We introduce the first forward difference operator acting on x 

=f{x+\) -f{x) 



as a discrete analogue to differentiation. Since 



£ 

;=n 



Af n (t)=f(x+i)-m, 



(a discrete analogue of the Fundamental Theorem of Integral Calculus) the first forward 
difference operator is a discrete inverse of the summation operator. This difference op- 
erator is used in the Rodrigues' formula for the Hahn polynomials, which can also be 
written in terms of the weight function w(.v;a, /?, N) with shifted parameters as was 
done for the Jacobi polynomials. 

Two limiting cases of the Hahn polynomials are of particular interest. In the 
first case, replace „y by Nx in the interval of orthogonality. This in effect places the 
support of the weight function on the ( A'+l) equally-spaced points { 0, 1/A', 2/A’, ... , 1 } 
in [0,1]. As N —* oo, the set of discrete mass points tends towards the full interval 
[0,1] ; we expect that this structure is reflected in the discrete orthogonality relation (22) 
becoming an integral orthogonality with respect to a continuous weight function on that 
interval. Rewriting the weight function in (23) as 



w(x ; a, /?, A r ) = 



r(.t4-l+oc) 

r(.v+l) 



r( A r — 1 +/?) 
r( A’-.v+i) 



and using a consequence of Stirling's Formula [Ref. 15: p. 257 j 
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we see that 



Hr+q) 

r(/) 



limf" 

oo 



r(a+ 1) r(/?+ 1) lim N' a ~ p w{Nx ; a, p, /V ) = (1-jc/, 

iV — ► oo 

which we recognize as the continuous beta distribution for the Jacobi polynomials nor- 
malized on [0,1]. 

This is indeed the case, and can be easily verified by a direct computation on the 
hypergeometric series definitions (21) of Q n {x ; a, /?, N ) and ( 1) of P { P' p \x) given in Sec- 
tion A.l, i.e., 



lim Q„( Nx ; a, p, N) = —^— p \l-2x). 

Thus the Hahn polynomials may be viewed as a discrete analogue and generalization of 
the Jacobi polynomials. [Ref. 14: p. 36] 

The second limiting case gives rise to an interesting class of polynomials which 
has applications in coding theory [Ref. 16]. 
a. Krawtchouk Polynomials 

For 0 < p < 1, let a = pt , p = {\—p) t in (21), then take the limit as t tends 
to infinity to obtain the Krawtchouk polynomials, i.e., 

lim Q n {x \pt ,{\ —p) t , N) = 

OO 

(25) ,F, ;}]- A> \p,N) 

[Ref. 14: p. 38]. A similar limiting argument applied to the weight function shows that 
the Krawtchouk polynomials are orthogonal with respect to a binomial distribution: 

(26) ^ K m {x ;p, N) K n {x ; p, K ) ^ ^ p x {\-p)‘ V ~ x = 0 , m¥=n. 

This class of polynomials possesses an inherent symmetry, the structure of 
which can be generalized to form other orthogonal polynomial classes. To this end, we 



56 



turn now to the characteristic of duality. For suitably defined functions p(x) and v(.v), two 
classes of orthogonal polynomials { p n (n(.x )) }, and { q n (v(x )) } are said to be dual if 

/>*(#*(«)) = 

That is, interchanging the roles of the degree n and the discrete variable x in one class 
produces the other. The Krawtchouk polynomials provide an example of an orthogonal 
class that is self-dual, i.e., 



K n (x-p, N ) = K x (n ; p, N). 



This is clear from the hypergeometric series definition (25) given above. It therefore may 
seem reasonable to suspect that there exists a class of orthogonal polynomials dual to 
the Hahn polynomials. Such a class does in fact exist, and it is this dual class which we 
next examine. [Ref. 17: p. 657] 

2. Dual Hahn Polynomials 

The dual Hahn polynomials, R n { ^.(jc) ; y, 5, iV), are defined 



(27) 



Rf ^ M > v, <5, A r ) — 3 F 2 



—n, — x , x+y+0+1 

y+F 



where /i(.r) = .v(x+y+<5+l), and n = 0, 1, ... , N. These polynomials satisfy the 
orthogonality relation 



(28) 



x 

'Y J R m { Kx) ; y, <5, N) R n ( /.(x) ; y, <5, A’) m{x ; y, (5, N ) = 0 , m # n 



x=0 



where 



, _ (y+<5+ 1), ((y+5+3)/2) x (y+ l) x (-N ) x (-if 

^ (Ox ((y+<5+ 1 )/2 )jc (<5-|- 1 )jc (y+<5+A'+2) JC 

_ . * (y+<5+ l) x (y+ l) x (— A~) x 2.v+y+5 + l 

( 1)^ (0T 1 ) x ( .V-hy+d+2) K y+t>+ 1 

Note that R n ( i(jc) ; y, 5, N) is a polynomial of degree n in the "variable" A(.v). The rea- 
son can be seen directly from the hypcrgeometric series 
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/?„( >.(x) ; y, <5, A* ) = 

/c=0 



(~ n )k U+y+^+Ofe 

(y+i),(-AO,/<! 



By writing the terms (— x)* (x+y+<5+l)* in product form, we obtain 



k - 1 k - 1 



J~| (—x+j) J~ [ (x+y+<5+/+ 1) 
y= o y=o 



which becomes 



A-i 

Y\ (~ x +j) (x+y+S+l+j). 
y=o 

Multiplying these factors as binomials, we have 



k - 1 

["] [ (-*) (x+y+5+l) +y'(-x) +y (x+y+<$+l) +j 2 ] 
y=o 



which simplifies to 



fc-i 

m -x (x+y+<5+l) +/(y+<5+l) +/ ]• 
y=o 

Taking only that part which depends on x yields 

/(x) = x(x+y+<5+l) 



as given above. [Ref. 18: p. 48] 

As discussed in the previous section, the discrete classes interact more naturally 
with difference operators than with differentiation. To accomodate the quadratic form 
of /(x), we introduce the divided difference operator 



Pf n { >.{*)) 



a a /(x)) 

a ;.(x) 



/„(;.(x+ 1 ))-/„( ;.(x)) 
;.(x+i)-;.(x) 



The dual interplay between the dual Hahn and Hahn polynomials follows easily 
from their hypcrgcomctric definitions: 
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= <?„(*; a, P, jV). 



R x { ;.(/;) ; a, p. 




—n, n+tx+p+\ 
a+1 ,-N 



This duality is also reflected in the recurrence relations and difTerence equations of the 
two classes. By interchanging the roles of n with x, a with y, and P with <5, the recur- 
rence relation for one class leads to or can be extracted from the difference equation for 
the other class. [Ref. 14: p. 37] 

3. Racah Polynomials 

Both the I lahn and dual Hahn polynomials can be unified as special cases of a 
single larger class, the standard notation of which is similar to that of the dual Hahn 
class. The Racah polynomials, R n { ^.(x) ; ot, p, y, (3), are defined by the hypergeomctric 
series 



(30) 



R n ( ?.(x ) ; a, p, y, <5) = 4 F 3 



~n, n-\-<x-\-p Tl, — x, .v+y+(5+l 
a+1, /?+(5 + l, y+1 



where i(ar) = ;c(x+y+<5+l), n = 0, 1, ... , N, and one of a+1, /?+( 5+1, or y+1 equals —N. 
(Physicists understand these objects via 6-j symbols.) 

The orthogonality relation for the Racah polynomials is 



,v 



(31) ^ R m ( ;.(.v) ; a, /?, y, 3) RJ A(x) ; a, /?, y, 3) w(x ; a, p, y, 3) = 0 , m ¥= n 



x=0 



where 



(y+<5+ IQ. ((y+<5+3)/2) JC (a+lQ ( P+3+ IQ (y+l) x 
(Ox ((y+<5+l)/2)jf (y+<5— a+0x (y — /?+0x (^+0x 
_ (y+<)+ Ox (fl+ 1 )x (/?+<)+ l) x ('/+ 0v 2.v+y+()+l 
( 1 Ij. (y+<5— a+ 1 ) x (y-p+l) x (<5+ 1 ) x y+3+ 1 

[Ref. 19: p. 24], Note that ; a, /?, y, c5) has a "well-poised" structure. This means that 
the pairwise sum of numerator and denominator parameters is constant, i.e.. 
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y+5+3 



(y+5+1) + (1) = y+<5+2 



+ 



y+5+1 



= y+<5+2 



(a+1) + (y+<5— a+1) = y+<5+2 
(/?+<5+l) + (y— /?+l) = y+<5+2 
(y+1) + (<5+l) = y+<5+2. 



The same could be said for the weight function of the dual Hahn polynomials. (To say 
that a hypergeometric series is well-poised means that for 

a Jf , a r+1 

h • v 



the parameters a and b satisfy the following relation 

tf]+l = = = Q r +\~bb r 



Wcll-poisedness is an important property of certain summable hypergeometric series.) 

From the hypergeometric series definition (30) of the Racah polynomials, we 
find that when the roles of * and n, a and y, and /? and <5 are all interchanged, the series 
is unchanged. Thus, like the Krawtchouk polynomials, the Racah polynomials are self- 
dual. 

To recover the Hahn polynomials as a limiting case of the Racah polynomials, 
let y+1 = —N and <5 -» oo. Thus by formulas (30) and (21), 



— n, n+a+P+l, —x, Jc+y+(5+ 1 




—n, n+a+P+l, —x 


a+1, /? +<5 + 1 , — A _ 


= 3 F 2 


a+1, -N 



and then by formulas (32) and (24), 



lim w(„v ; a, /?, y, «5) = 

6— ►oo 



(0,v 

(/?+D,V 



w(x;a, P, A'). 



Likewise, by letting a+1 = —N and /? -> oo, we obtain the dual Hahn polynomials. 

Figures 5 and 6 together with Tables 2 and 3 provide a hierarchy of classes dis- 
cussed in this chapter. Tables 4-15 summarize information about selected classes in this 
hierarchy [Refs. 14, 12, 20, 21 ,5, 13, 9 , 22]. 
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Jacobi 



Gegenbauer 




Chebyshev Chebyshev Legendre 

First Kind Second Kind 



Figure 6. Hierarchy of Special Cases 
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Table 2. TABLE OF CLASSES 



Class 


Symbol 


I Iypergeometric Series 


Racah 




r V—n, n + a + fi + 1 , — .r. „v + y + <5 + 1 ,1 
4/ ' 3 L a+l,/? + <5+l,y+l ’ ‘J 


Dual Hahn 


R n ( A(x) ; y, 5, N ) 


p V-n, ~ x > x + a + /? + 1 . i ”1 
3 /2 L V + 1. -iV (1 J 




Hahn 


Q n ( x ; «> P, N) 


r- f— x, — n, n + a + P + 1 . 

+ a + 1, -iV ’ J _ 




Jacobi 


Pn'%) 


( a + 0/t r T—a/, n + a + /? + 1 l-.r 1 

n! 2 '>L a +! ’ 2 J 


Krawtchouk 


K n {x\p, N) 


p r _ _ • — ”i 

2^1 l _ jY > p J 


Meixner 


M n (x ; P, c ) 




Laguerre 


r(*), s 

Ln M 


(“ + 0/1 r. r -n 1 


Charlier 


QC* ; a) 


V-n, —x . L "I 

2 y o[_ _ - a J 


Hermite 




(2r) n 2 F 0 ^ -/,/2 ’ ( + O/ 2 . _ _L 
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Table 3. TABLE OF CLASSES ( 


CONTINUED) 


Gegenbauer 


c n M 


(2^)« r V —n, n + 2/ 1 — jc "1 

n\ 2/, L /+1/2 ’ 2 J 


Chebyshev, 
First Kind 


U.x) 


1 1 

1 ^ 

£ 

1 

i i 

t-C 

CM 


Chebyshev, 
Second Kind 


U n (x) 


f i i \ p r tl % H 4" 1 1 X j 

(«+l)2^i[ 3/2 ’ 2 J 


Legendre 


P n (x) 


F r~«, n+ 1 . 1 - .r “1 

jF ‘L 1 ' 2 J 
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Table 4. RACAH POLYNOMIALS 

Symbol: R„( AO) ; a, (l, y, (5) 



Interval: For N a positive integer, x = 0, 1, ... , N. 

Weight: 

(y+<5+l) x ((y+<5+3)/2) x (a+l) x ( /?+i5+l) x (y+l) x 



wO ; a. P, V, <5) = ' 



(l)* ((y+o+l)/2) x (y+<5— a+l) x (y—p+l) x (o+l) x 



Norm: 



N 



2] [ R n( A M ; a, P, y, S) f w(x ; a, p, y, (5) 

n\ {n+o.+p+ 1)„ ( /?+!)„ (a— <5+l)„ (a+/?— y+l)„ 



jc=0 



= M 



where if 
a+1 = —N, 



(a+/?+2) 2 „ (a+l)„ ( P+5+ 1)„ (y+l)„ 



, . _ (y+<5+2) / y(— P) N 

' " (y-/?+l) lV (<5+lO ’ 

or if 

P+S+l = -N, 

_ (y+<5+2), v (<5— a),v 

(y+(5— a+l) /V (<5+I) iV 

or if 

y+1 = —N, 

. . _ (~<5),y (° t +/?+2), v 

(a— t)+l) (V ( /I+l)y 
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Table 5. RACAH POLYNOMIALS (CONTINUED) 



I Iypergeometric Series: 



R n ( /(.y) ; a, /?, y, <5)= 4 / 7 3 



P — u, n+a+P+l, — x, .Y+y+(5+l 
oc+1, p+3+\, y+1 




where l(.v) = x(x+y+()+l), n = 0, 1, ... , N, 
and one of a+1, /?+(5+l, or y+1 equals — N. 



Recurrence Relation: 



/(.v) R n { i(jc) ; a, /?, y, 5) = a{n) [ y.(.v) ; a, /?, y, ;.(x) ; a, /?, y, <5) ] 

-c(n) [ R n ( ;.(jr) ; a, p, y, 3)-R n _P />.(*) ; a, p, y, (5) ] 



where 



{n+a+P+l)(n+a+l)(n+P+3+l) (rt+y+1) 
(2n+a+P+l) (2n+a+p+2) 
n ( n+P ) («+a+/?— y) (n+a— 3) 

(2«+a+/?+l) (2«+a+/?) 



Difference Equation: 

;.(«) /?,,( ;.(«) ; a, P, y, 3) = A(x) [ /?*+,( ).(n ) ; a, 

-av)[^( ;.(«); a, 



P, y, <5)-Z? x ( yl(«) ; a, P, y, <5) ] 
P, y, <5)-/?^_,( A(rt) ; a, p, y, (5) ] 



where 



4 *) = 



c(*) = 



(.Y+y+d+l) (.v+y+1) (^-+<3+/?+ 1 ) (.v+a+1) 
(2.r+y+(5+l) (2.Y+y+<5+2) 

.y (.y+< 5) (x+y+h— a) (x+y—£) 

(2.v+y+(5+l) (2.r+y+(5) 
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Table 6. DUAL HAIIN POLYNOMIALS 

Symbol: R„( A(jtr) ; y, 3, N) 



Interval: For N a positive integer, .v = 0, 1, ... , /V. 

Weight: 

, . y . (y+(5+l) x ((y+<5+3)/2) x (y+l) -t (— N ) r (— 1)* 

(1)jc ((y+c)+ l)/2) x (<5+l) x (y+<5+A , +2) JC 



Norm: 



X t ; -W ; <5 ’ N ) ] 2 w; (- r ; y* ^ ^ ) = ( /z „ y ) ( A / v-« <5 ) 1 ' 

jc=0 L J 

Hypergeometric Scries: 

R„( >.(x) -,y,i,N) = 3 f 2 [- n ' 1 ; l] 



where y.(.v) = x (jc+y+5+1) 
Recurrence Relation: 



-;.(*) R„( A(x) ; y, <5, N) = D{n) [ R n+l ( ?.(x) ; y, (5, N)-R n ( ?.(x) ; y, (5. N) ] 

—D(n) [ R n ( ; y, <5, N A(x) ; y, <5, iV ) ] 



where 



Z?(«) = ( A — /?) (y+ 1 +//) 
D{n) = n{ N+ l+<5— n) 



Difference Equation: 

-nR n { ).{x) ; y, 5, N) = b(x) [ R n ( X(x+l) ; y, <5, N)-R n ( ).(x ) ; y, <5, A')] 

-d(x) [ R n ( A{x) ; y, 3, N )-R n ( y.(.r-l) ; y, 3, N ) ] 

where 



_ (.v+y+b+l) (.r+y+1) ( ,V-.v) 
X (2w+y+<)+ 1 ) (2.c+y+<)+2) 

_ .r fv+(3) (.v+y+(>+.V+ 1 ) 
(2.v+y+<3) (2.v+y+«)+l) 
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Table 7. HAHN POLYNOMIALS 



Symbol: Q„(x-,a, P, N) 




Interval: For N a positive integer, x = 0, 1, .. 


.,/V. 


Weight: 

e. »)-(*?) ( n n x -T, 


i (A+O/v (a+ 1 ) x (— iY) x 
(i),v (i) A-x-P) x 


Norm: 




oo 

^ [ Q„{x ; a, P, N ) f w(x ; a, /?, /V ) 

x=0 

a+ZM-1 


(N+a+p+l^ (N+a+p+l+n^ 


2/7+a+/?+l 


0 


r(a+i)r(a+/?+i) n«+/?+i)r(«+i) 


r(/?+i) 


r(«+a+l) r(/?+a+/?+l) 


Hypergeometric Series: For a, P > — 1 and 
n = 0, 1, ... , N, 




&(*;«. /». 'V)- 3 F J [-"-"+ + “+i+ l '- Jr ; 1 ]. 
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Table 8. HAHN POLYNOMIALS (CONTINUED) 



Recurrence Relation: 

-xQ n 0;a, P, N ) = &(«)[ Q, i+1 (jt; a, /?, A r ) - ; a, P, /V)] 

-4«) [ Qn(* \«.P,N)~ Q n - ,(* ; a, P, iV ) ] 

where 

(n+a+/?+l) (tf+a+1) ( A— «) 



*(«) = 
d{n) = 



(2n+cc+p+ 1) (2«+a+/?+2) 
n(n+p) (n+a+P+N+l) 

( 2n+a+P ) (2//+a+/?+l) 



Difference Equation: 

-;.(«) &,(* ;a, P, N) = B(x) [ &,(*+ 1 ; a, p, N)- Q n (x ;«,/?, iV ) ] 

-Z)(.v) [ &(* ; a, P, N ) - &(x-l ; a, /?, iV ) ] 

where 

Z?(.t) = ( N—x) (a+l+.v) 

D(x) = x ( N + 1 +p—x) 

/(//) = h (/7+a+/?+l) 



Rodrigues' Formula: 

(T)C v ,v-^)('D aU; “' AlV) 



_ /«+/?\ * n fx+a\ ( N—x+P+n\ 

- v « / A W w v /*+« y 
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Table 9. JACOBI POLYNOMIALS 

Symbol: I*' p \x) 



Interval: [ — 1,1] 

Weight: (1 — x) % (1 + x) 1 * 

Standardization: I*' P> (1) = ( n + a ) = — + 

\ n J n\ 

Norm: 




[P ( n' ^V)] 2 (1 - xf (i + xf dx = 



2 a+/?+1 r(«+«+l) rjn+fl+l) 
(2n+a+fi+\) n\ r(tf+oc-l-/?+l) 



Hypergeometric Series: P ( n '' p \x) = 2 ^ n ’ "+ a +^ +1 \ yy J 

Recurrence Relation: 2 (n+l) (n+a+fi+l) (2«+a +/?) P^i p \x) 

= (2/?+a+/?+l) [ (2 n+a+p) (2«+a+/?+2) .v+a 2 -/? 2 ] P l *'%) 

-2 («+a) (2/i+a+/?+2) P { *lf\x) 



Differential Equation: 

(l-jr 2 )/' +[/? — « — (oc+/?+2) Jf] /+« (/7+a+/?+ 1) >> = 0 , y = P^' P \x) 
Rodrigues' Formula: 

2” *! P ( :~ %) = (- 1 ) n ( 1 ( 1 +-yf /? -pr [( 1 - *) a+ "( 1 + 

dx 

Generating Function: 

CO , o 

y P„’ p \x) t n = — == ( \—t+yJ l - 2xt + r 2 V* (l+r+^1 - 2.vr + r 2 

n=0 J(l-2xt+t 2 ) V 7 V 7 

Explicit Expression: /f "V) = 2'" £ ("J") (jc- 1)-* (.v+1)* 
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Table 10. LAGUERRE POLYNOMIALS 

Symbol: L<%v) 

Interval: [0, oo) 

Weight: .v* e ~* , a > — 1 



r_ i v 

Standardization: L { *\ 1) = — — — 



n\ 



.r+ 



Norm: 



j n n. 



Hypergeoinetric Series: 



h’V)-(Tk.U'H 



Recurrence Relation: 



(rt+1) L*+,(jc) = [ ( 2n+a+ \ )-x ] L^’(x) - ( n+a ) L^\{x) 



Mr 






DilTcrential Equation: x y" +( a+ 1— x)y'+ny = 0 , y = L { *\ x) 
Rodrigues' Formula: 

r (a)/ \ 1 dr >t+a ~ x 1 

L\\x)= , _, x —ir[x e ] 

n\ x e ax 



Generating Function: 



OO 

^ exp( 7 f r ) 

n=0 



Explicit Expression: 



OO 

^w = E(-o*(S)^,‘ 



0 
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Table 11. HERM1TE POLYNOMIALS 

Symbol: //„(*) 

Interval: (— 00 , 00 ) 

Weight: e ~* 2 

Standardization: // n (l) = T x n + ... 

Norm: 




—OO 



llypergeometric Series: 




2 



Recurrence Relation: 



H n+] (x) = 2.x II n (x) - In //„_,(*) 



Differential Equation: 



y" ~ 2 xy' + 2 ny = 0 , y = II n (x) 



Rodrigues' Formula: 




Generating Function: 




Explicit Expression: 




D/2] 
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Table 12. GEGENBAUER (ULTRASPHER1CAL) POLYNOMIALS 

Symbol: C\x) {or P?\x)) , A > —1/2 

Interval: [—1,1] 

Weight: (l-.v 2 )^ 2 
Standardization: 



C<%) 

lim — — — 

Cf(l) 



= T„(x) , n = 0, 1, 2, ... 



Norm: 



J 



' \C {A \x)V (l - x 2 )' 1-1/2 dx = — 

L n 1 > («+;.)«! n«+2;.) 



Hypergeometric Scries: 



- r —n , n+ 2X 1— .V "1 

'»L 2+1/2 ; 2 J 



W rt ; 2M I ;.+i/2 ’ 2 



Recurrence Relation: 



(« + 1) C$,M = 2 («+^) x C*V) - («+2i-l) C-! W 



We 



-+■) 



DilTerential Equation: 

(1-JC 2 )/-(2A+1) -X /+/1 («+2;.)>’ = o , y = C ( n \x) 
Rodrigues' Formula: C 0 w (.v) = 1 , C ( /\x) = 2).x, 

2 n n\ U+l/2)„ (l-x 2 y-' 12 Cf (*) = (-1)” (2 >X-£r[ {l-x 2 T +A ~ U2 ] 



Generating Function: 



o o 

E 

n=0 



Cf(.v) t n = (l-2.vr+/ 2 )- 
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Table 13. CHEBYSHEV POLYNOMIALS OF THE FIRST KIND 



Symbol: T n (x) 

Interval: [—1,1] 

Weight: (1 — „v 2 )" 1/2 
Standardization: 7 n (l)=l 
Norm: 

V,wf( i— 2 r ,,2 *={f';;fo 

*^-i 



Hypergeometric Series: 




] 



Recurrence Relation: T„ +l (x) = 2xT n (x) — T„.,(x) 

DifTerential Equation: (1 — .t J ) — xy' + n 2 y = 0 , y = T n (x) 
Rodrigues' Formula: 



T( \= - ilil * I 1 r/i 2\n— 1 /2"l 

nX 2 n+] r(« + 1/2) dx n ^ ^ ] 



Generating Function: 

oo 

n =0 



1 —Xt 
1—2 xi+i 2 



— 1 < jc < 1 , | / | < 1 



Explicit Expression: 



T n (x) = cos (n9) with x = cos 0 , 0 <, 0 <, n 
[«/ 2] 

r "W*f 



-n» 



*=0 



A!(/t - 2A)! 
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Table 14. CHEBYSHEV POLYNOMIALS OF THE SECOND KIND 

Symbol: U n (x) 



Interval: [—1,1] 

Weight: (l-.r 2 )'/ 2 
Standardization: t/„(l) = «+l 
Norm: 



J 



c^w] 2 (i - x 2 y l 2 dx = —■ 



Hypergeomctric Series: 



C/„( J r) = (/i+l) 2 fi[ rt j/2 + 1 ; ±JL 



Recurrence Relation: U nM (x) = 2xU n (x) — i/„_,(.v) 
Differential Equation: (1 — x 2 )y" — 3 xy' + n {n+2)y = 0 , y 
Rodrigues' Formula: 

(-l)' 1 (/7+ \)Jn t f 



U n {x) = 



(l-.v 2 ) 1/2 2' ,+1 r(/i+3/2) dx 



-[(1 



Generating Function: 

oo 

Z = 1 



n= 0 



1—2 xt+t 



, -1 <* < 1 , | 



Explicit Expression: 



sin(rt + 1)0 

U n i x ) = TT, with * = cos 0,0: 



sin 0 



D/2] 

v.w = X! < -d‘ 



in — A)! 



k = 0 



A! (// - 2A)! 



i2x) n 



= U n (x) 

■x 2 y +] ' 2 ] 

t\ < 1 

S 9 < 71 
-2 k 
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Table 15. LEGENDRE (SPI 1ER1CAL) POLYNOMIALS 

Symbol: P„{x) 



Interval: [—1,1] 

Weight: 1 

Standardization: P n (l)= 1 



Norm: 



fl 



[/\(x)]Vx = 



2 

2/i+ 1 



Hypergeomctric Series: 



PM 




—n, n + 1 

1 



i-x 

2 



] 



Recurrence Relation: {n+ 1) P n +M = (2/i+ 1) x P n {x) — n P„.,(x) 
Differential Equation: (1 — x 2 )y" — 2xy‘ + n (/i+ 1) j/ = 0 , y = P„(x) 

Rodrigues' Formula: P n (x) = ^ [ ( 1 -x 2 )"] 

Generating Function: 

OO 

^ /^(x) t n = (l— 2 xH -/ 2 ) _i/2 , — 1 < x < 1 , | / | < l 

n~ 0 

Explicit Expression: P„(x) = -jr g ( -1)* (£) ( 2n ~ 2k )^~ ik 
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V. APPLICATIONS 



Orthogonal polynomials and special functions in general have been studied exten- 
sively in mathematics and other fields since the eighteenth century. Presented below are 
a few of the traditional applications of selected classes of orthogonal polynomials. Our 
first few applications come from numerical analysis. 

A. ECONOMIZATION OF POWER SERIES 

Economization of power series is a technique used to reduce the degree of a 
polynomial approximation to a given function. 

The maximum norm (or L°° — norm) for a continuous function on a compact inter- 
val \_a,b ] is defined as 



This norm is not induced by an inner product, but nevertheless has many uses in applied 
mathematics, including numerical analysis. 

The minimax property of the Chcbyshev polynomials states that of all n* degree 
monic polynomials (i.e., leading coefficient 1), 2'~ n T n (x) has the smallest maximum norm 
on [ — 1,1] [Ref. 23: p. 106]. The justification for this statement is deferred until Section 
C. Hence the best approximation in the maximum norm to the function x" on [ — 1,1] 
by a function of lower degree is f(x) = x n — 2'~ n T n (x). So, given a function and a 
polynomial approximation to that function (e.g., from a Taylor series expansion), suc- 
cessively replace the highest powers x" with f(x) to obtain a polynomial approximation 
of lower degree [Ref. 23: p. 125j. 

Example 1: Let f(x) = sinx. The Maclaurin series for this function is 



and this series is convergent for x e (R. If truncated after the x 5 term, the polynomial 
approximation for the function is 



ll/IL = max I/M I • 

xs [a, b] 




sin x ~ x — 




77 



with a maximum error of 0.0002 for are [—1,1]. We use the fifth degree Chebvshcv 
polynomial T s (x) = 16M — 20„v 3 + 5x to obtain 



,r 5 ^ 5.v 3 _ 5x 
5! * 4 SO “ 1920 



with an error not exceeding 



max 1 r s (.r) | _ i 

1920 “ 1920 



0.00052 



in [ — 1,1]. Thus the approximation 



sin x « x — 




+ 



< 3 
5x 

4S0 



5x \ = 1915 _ _75_ 3 

1920 ) 1920 X 4S0 X 



has an error whose magnitude in [ — 1,1] does not exceed 

0.00052 + 0.0002 = 0.00072. 

Compare this with the maximum error of 0.00833 for the Maclaurin series which is 
truncated after the x 3 term. For a cubic polynomial approximation of sin x, the "econ- 
omized" polynomial has a maximum error that is significantly smaller (less than one 
tenth) than that of the truncated Maclaurin series. 

The next three applications illustrate the usefulness of the zeros of orthogonal 
polynomials. The first two come from numerical analysis, the third from a problem in 
electrostatics. We begin with a preliminary discussion of a fundamental technique from 
numerical analysis. 

B. POLYNOMIAL INTERPOLATION 

Polynomial interpolation is a method of approximating a given function with a 
polynomial that matches ( interpolates ) the function at specified points (called nodes or 
abscissae) .t,, ... , x n (Ref. 24: p. 497], Given a function /(.v) and n distinct nodes in a 
compact interval [a, A], there is a unique polynomial of degree (n — 1) that passes 
through the points (.v, , f(xj), l <, i < n. 

For each /= 1,2, ... , n , define a polynomial of degree (n — 1) by 
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- x ~ x k 

*l-*k 



k±l 



Clearly, for each k = 1,2, ... , n, n itn (x k ) = 0 if /< # / and 7r (/1 (.v,) = 1, that is, 7r lV ,(.tJ = <5,*. 

n 

(We may thus equivalently express ^(.r) = l n (x)j{{x—x) l n '(x)), where l n (x) = f [ (w— x k ) is 

A= 1 

the unique monic polynomial of degree n that vanishes simply at each node x k , 
k= 1,2, ) 

The Lagrange interpolating polynomial is now given by 



n 

L f ni x ) = T/Ct) Xifii*)- 
/= 1 

Clearly, L{(x k ) =f(x k ) for k = 1,2, ... , n, and uniqueness of the degree ( n — 1) interpolat- 
ing polynomial is guaranteed by the Fundamental Theorem of Algebra. 

C. OPTIMAL NODES 

In this section we address the issue of estimating the maximum size of the interpo- 
lating error || f ~ L f n \\ x . Assuming f is suitably differentiable in \_a,b~\, it can be shown 
[Ref 25: p. 188] that there exists a value K n (which depends on /'"’in [a, 6]) such that for 
any x in [o,Z>], 



I /(.*)- £.{(*) I <-§- KM I, 

n 

where l n (x) = fl ( x ~ -F)- (Note again that l n (x) is a monic polynomial of degree n which 
1=1 , 

vanishes at the nodes.) This implies that \\f— Li |L is minimized by making the optimal 
choice of nodes jq, ... , x„ in [a,b] which minimizes || l n U^,. Surprisingly perhaps, this hap- 
pens precisely at the zeros of the Chebyshev polynomials T„(x), scaled to the interval 
We sketch the reasons below. 

For simplicity, we take our interval of interest to be [ — 1,1] instead of [u,Z>] without 
loss of generality. The transformation 



x = 




or equivalently, 
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'= 2 (i£r)-'- 



is a one-to-one continuous mapping between the intervals [ — 1,1] and [a,Z>]. 

Now let .v, be such that T„(x t ) = 0, / = 1, 2, ... , n (see Chapter I, Section A. 3), i.e., let 
l n {x) = ~^—p T n (x). From the recurrence relation (Equation (4) in Chapter I), r„(jt) has a 
leading coefficient of 2 n_l ; hence this l„(x) is monic. Moreover, from the definition 
T„{ cos (?) = cos nd, it follows that 

a) t.w-t-o'-jb- 

where y, =y tJ , = , / = 0, 1 , ... , n. Now suppose m n (x) is another monic polynomial 

of degree n such that 

(2) ii io» < ii u loo - "pr - 

Combining (1) and (2) we see that we must have for / = 0, 1 , ... , n, 

— if / is even, and 

I* 

— V if i is odd. 

yi-i 

jL 

Thus, the polynomial p n (x) = m n {x) — l n (x ) has degree at most (n — 1) (since it is a differ- 
ence of two monic polynomials of degree n) with at least n zeros, one in each interval 
(y, , 0 < i <n — 1, by (3). This contradicts the Fundamental Theorem of Algebra, 

and so no such polynomial m n (x) satisfying (2) exists. Hence the choice of nodes jc, de- 
termined by l n (x ) = - * 7Xv) minimizes || /„ = max I l n {x) \ over all possible monic 

polynomials l n (x), and therefore over all possible choices of interpolating nodes x r 

It should be emphasized that this Chebyshev interpolation allows an a priori error 
bound for all .t, but is not always best possible for every .v using other interpolation 
schemes. (For example, if equally spaced nodes „v„ ... , x„ are used, then trivially, the er- 
ror f(x ,) — L/ n (x,) = 0, even if .v = .v, is not a Chebyshev zero.) Chebyshev zeros arc opti- 
mal when one has freedom in the choice of nodes. For details regarding the practicality 
of Chebyshev interpolation and some asymptotic results, see [Ref. 24]. 



( 3 ) 



m n(y<) < Uyi) 

m n(yi) > 



80 



D. GAUSSIAN QUADRATURE 

Quadrature formulas are used in that area of numerical analysis concerned with the 
approximate integration of a function f{x) against a weight function vv(jc) > 0 on an in- 
terval ( a,b ), when the explicit evaluation is intractable. 

An inter polatory quadrature is such a rule that uses interpolating polynomials, such 
as Lagrange polynomials: 



[ fix) w(x) dx » f L f n {x) u--(.r) dx, 

** a ^ a 

n 

where L f n {x) = Y,fi x i) n ,(x), as described above. This can be rewritten as 

f=i 

(1) J f{x) wW dx « W„ 



where the weights w, are given by 



(2) W[ = f 7t i n (x) w{x) dx. 

^ a 

Thus, interpolatory quadrature is basically a weighted sum of the function values f{x) 
at the nodes x, , /= 1,2,... ,n, as are numerical integration recipes such as Simpson's 
Rule and the Trapezoidal Rule. 

For specified nodes .t,, ... , x n e the n weights w u ... , w n computed in (2) for the 

quadrature (1) will be exact for polynomials / e Pn-faf], but we can do better. [Ref. 
25: p. 236] 

In Gaussian quadrature, we ask for the location of the n nodes jr„ ... , x„ as well as 
the n weights w 1( ... , w n in order for the quadrature rule (1) to be exact for polynomials 
f e P 2r ,.i[a,b]. At first glance, this seems to be an extremely complicated computational 
problem, but the solution falls out simply when the theory of orthogonal polynomials 
is applied. 

Let { p„{x) be the class of polynomials orthogonal with respect to the weight 
function vv(.v) on [a, 6], say, by Gram-Schmidt, and let iff. x) be the Lagrange polynomial 
that interpolates / e P 2n _ ,[«,/?] at the zeros of pfx), so that f(x) — If(x) e P 2n _S_a,b~\ van- 
ishes at the zeros also. Since the zeros of p„{x) are real, simple, and lie in («,Z>), we have 
the property that 
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fix) - I&x) = Pnix) ? n -lW. 



where e P n _ ,[«,£>]. l ienee, by orthogonality, 

f f{x) w(.r) dx - [ lf n {x) w(w) dx = [ p n (x) q n _fx) w(x) dx = 0. 

J a J a J a 

We have thus shown that the resulting n-point Gaussian quadrature rule is exact for 
/ e P to _|[a,6]. Moreover, by a theorem of Sticltjes, if f(x) is continuous on a finite inter- 
val [<?,£>], then 



lim 

n ~*oo 



'b rb 

L f n {x) vv(x) dx = I fix) w(.v) dx. 



By applying the Christoffcl-Darboux formula (Equation (9) of Chapter III with 
notation from Chapter III, Section B) and using (1), it is possible to derive an alternate 
expression for the weights (2): 



^n + 1 

w- = — 

K Pn( x i) P n +li x i) 

which are referred to as the Christ off el numbers. [Ref. 8 ] 

Thus Gauss-Jacobi, Gauss-Chebyshev, and Gauss-Legendre are the names given to 
Gaussian quadratures involving the weight functions and orthogonal polynomials from 
the Jacobi, Chcbyshev, and Legendre classes, respectively. 

E. ELECTROSTATICS 

The zeros of the Jacobi polynomials play an interesting role in a problem of Sticltjes 
concerning electrostatic equilibrium. In this problem, fix "masses" of positive charge a 
and P at the points .v = 1 and x = — 1, respectively. Then place n point masses of positive 
unit charge in the interval ( — 1,1) so that they arc free to move. These interior masses 
arc now subject to a "logarithmic potential", that is, a repelling force that is proportional 
to the logarithm of the distance separating them. The problem is to determine the dis- 
tribution of the point masses x , , /= 1,2, ...,«, when the system is in equilibrium. 
Mathematically, this is equivalent to maximizing the force function 
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/— 1 /=1 l<i<j<n 



The logarithmic terms give the restrictions j:,#— 1, x n ¥= 1, and x,#x ; - for ij=j. This 
otherwise continuous function gives the equilibrium points by setting 



( 4 ) 



t- = 0 , /= 1, 2, ... , n 

ox t 



and solving for x u ... , x„. To solve this system of n nonlinear equations in n unknowns, 
Stieltjes introduced the polynomial 



ri 

Prk x ) = P] ( x ~ x i) 
1=1 



and reduced (4) to 



1 P 

2 Pn'( x i ) {+x i 



1 Xi 



= 0 , 



1,2, ... , n 



which becomes 

(l-x 2 >„"W + [ 2p - 2a - (2a + 2/?) x t ] p n '(X [ ) = 0 

for / = 1,2, ... , «. 

Since the polynomial 

(l ~ x 2 ) Pn\ x ) + [ 2/? - 2a - (2a + 2/?) x ] p n ’(x) 

is of degree at most n and vanishes at x = x„ / = 1,2, ... , n, it can be set equal to a scalar 
multiple ).p„{x ) which also vanishes at these points. Hence 

( 1 “ * 2 ) P n "( x ) + [ 2p - 2a - (2a + 2/?) x ] p n '{x) - ).p n {x) = 0. 

Attempting a power series solution to this second-order differential equation leads to the 
observation that polynomial (i.e., terminating) solutions exist if and only if 

A — — ii(n + 2 a + 2/? — 1). 
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Rearranging terms yields 



(1 - * 2 ) p~{x) + [ (2p - 1) - (2<x - 1) - ((2a - 1) + (2/? - 1) + 2) * ] p n '{x) 
+n(n + (2a - 1) + (2/? - 1) + 1) p n {x) = 0 



which is the differential equation for P^~ u2p ~ l) (x). Thus the equilibrium positions of the 
unit charges occur at the zeros of this Jacobi polynomial. [Refs. 8, 13: p. 140] 

The zeros of the Laguerre and Hermite polynomials can be developed as solutions 
to similar electrostatic equilibrium problems. 

F. SPHERICAL HARMONICS 

We investigate another application of orthogonal polynomials to problems in 
mathematical physics. In Cartesian coordinates (x,y,z), the Laplacian operator of a 
function u(x,y, z) is defined as 



^2 ^2 
c a , o a . d u 

V u = — r 4 r + 



dx 



2 



dy 



dz 



2 



When converted to spherical coordinates (r, 0, 0), this operator acting on a function 
u(r, 6, 0) becomes 



( 5 ) 




(■ r 2 u r )r + — TW ( sin 0 “e)e + ~ 

oil! U ci 



sin 2 0 



A e e 



where subscripting with r, 0, or 0 denotes partial differentiation with respect to that 
variable. A function u is said to be harmonic in a region D if in that region it satisfies 
Laplace's equation: 

V 2 u = 0. 



In particular, if the boundary of D is the unit sphere centered at the origin, then using 
the method of separation of variables, the Legendre polynomials arise naturally as part 
of the solution. 

In separation of variables, we assume that the solution will be of the form 
(6) u(r, 0 , 0) =f{r) g{0) /i(0). 

Substituting (6) into (5) and dividing through by u, we obtain 
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( 7 ) 



1 






h sin 0 



2 n h<p<p 9 . 



The first term, depending only on r, must reduce to a constant which we write as 
v(v+l). Substituting v(v+l) into (7) and multiplying through by sin 2 #, we have 



v(v+ 1) sin 2 # + ( sin # g 0 ) 0 + y h 4 , <p = 0. 



Now we see that the third term, depending only on <f>, must also reduce to a constant; 
call it —in 2 . Substituting —m 2 for this term and simplifying, we obtain 



( 8 ) 



sin 2 #g0 0+ sin# cos # g 0 + [ v(v+l) sin 2 # — nt 2 ] g = 0. 



By the change of variable x = cos #, (8) becomes via the Chain Rule 



/, 2\ d S 

l 1 "* ) TT 



— 2x-f- + 
dx L dx 



v(v+l) — 



m 



\-x l 



g = 0. 



For problems that are radially symmetric, the 4> dependence can be removed by setting 
m = 0, leaving 

,2 , 

(9) (l-w 2 ) — § 2x-^- + v(v+l)g = 0. 

When v is a positive integer n, (9) is recognized as the dilferential equation for the 
Legendre polynomials. [Ref. 26: pp. 210-213] 

An alternate approach to this problem (by the method of images) uses Green's 
functions. As motivation, consider the Dirichlet problem for the unit circle in the plane, 
which involves finding a harmonic function u(r, #) in the unit disk that takes on pre- 
scribed function values /(#) on the boundary r= 1. The solution is given by 



u(r, #) 



1 f 2n 

' 2* J o J 



f(t)P r (0-0') d0', 



where 



/)(>') = 



1 —r 



I — 2r cos y+r 1 
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is the so-called Poisson kernel for this problem. 

Similarly, the solution to the Dirichlet problem for ihc unit sphere in 1R 3 



W 2 u = 0 

u =f{0, 0) on r = 1 



can be expressed as 



u(r, 6, 0) = P7(0\ 4>') PM sin O' dO' df 



where 

( 10 ) 



PM = 



1 -r 



[ 1—2 r cos y+r 2 ] 3 ' 2 



and 



cos y = cos 6 cos O' + sin 0 sin O' cos(0 — 0')- 
Writing (10) in terms of simpler functions 



PM = 



Vl-2r 



cos y+r 



+ 2r ~ 

2 C7A* 



- V "l—2r 



cos y+r 



we note the appearance of the generating function for the Legendre polynomials, with 
x = cos y. Therefore the Legendre polynomials are again part of the solution. (Ref. 27: 
pp. 87-89] 

Laplace's equation can also be solved in a higher dimensional setting. Let 
x = (x H ... , x p ) e RL Note then that we may write x = ri, where 
r = || x || = fx 2 { + ... + x 2 p and ^ = — = (<;„ ... , l; p ) e R' is a unit vector. A polynomial 
h n (x) is said to be homogeneous of degree n if h n (/x) = A n h n (x). 

The Laplacian operator acting on a function u(x) in p dimensions is defined as 






1=1 



dji 

dxf 
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As before, a harmonic function is one which satisfies Laplace’s equation = 0. We 
now seek homogeneous harmonic polynomials of degree n in x e IR'. It can be shown 
that there are exactly 



linearly independent such solutions, and they can be characterized by Gegenbauer (or 
ultraspherical) polynomials C^tJ, k= 1,2, , N. The general solution is given by 



with q* suitably chosen unit vectors. Note that if p = 3, then these reduce to the 
Legendre polynomials found earlier. [Ref. 6: pp. 168-183] 

G. GENETICS MODELING 

Karlin and McGregor gave an interesting application of the dual Hahn polynomials 
to a model in genetics. In this continuous time Markoff chain model, the dual Hahn 
polynomials R n ( i(x) ; y, <5, N) arise in the transition probability function for the process. 

The setting for the model assumes N gametes of type a or A and gives a random 
fertilization scheme. The population of either type of gamete is affected by both the 
fertilization process and a mutation process whereby a gamete resulting from a mating 
can mutate into the other type. 

By considering the conditional probabilities for an increase in population size of 
both gamete types, a stochastic process is defined which is a classical birth and death 
process. The transition probability function for this last process is then cast in terms 
of the dual Hahn polynomials. The interested reader is referred to [Ref. 14] for the de- 




fyj(x) = K( r Q = r" S„(c), 



where the spherical harmonic 




tails. 
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VI. BASIC EXTENSIONS 



A. BASIC HYPERGEOMETRIC SERIES 

In this chapter, we extend the structure and some of the results in Chapter IV to a 
more general level. This extension is accomplished by introducing a new parameter 
called the base to the hypergeometric series. The base q was used by I Ieine in a series 



where c ± 0, — 1, —2, ... [Refs. 28,29]. This scries converges absolutely for | x \ < 1 when 
I q | < 1 by the Ratio test. Since 



we see that the series in (1) tends termwise to the ordinary hypergeometric series as 
q->\. Thus Heine's series is called the basic hypergeometric series or the q- 
hypergeometric series. [Ref. 12: p. 3] 

The q-shifted factorial is the basic extension of the shifted factorial introduced in 
Chapter II and is defined 




(l-<7) (l-? c ) (1-?) (l-V) (1-<7 C ) (l-<7 C+1 ) 



( 2 ) 



1 a 

\-q 



( 3 ) 




The ordinary shifted factorial is recovered by applying (2) and (3) in the limit 




We also define 



oo 



( 4 ) 
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a form vvc will see in later results. Since the infinite product diverges when both a¥= 0 
and ] q | > 1, we will assume | q | < 1 whenever (a ; q)^ appears unless otherwise stated. 
[Ref. 12: p. 3] 

Generalizing the basic hypergeometric series above, we define the ,4> s basic 
hypergeometric series (or r <p s series) 



(5) A 



a \ j ••• i a r 



b x ,...,b s ' q ' X J ^ {q\q) n {b ] \q) n ...{b s \q) n 



on 

=s 



(u, ; q) n ... (a r ; q) n 



(<-o\Q) 



1+w .r" 



//A 



where ^ J = n(n — l)/2 and q # 0 when r>s+l. We require that the parameters 
b u ... , b s be such that the denominator factors in each term of the series arc nonzero. 
Since 



(q m ; q)n = 0 - n = m + \,m + 2, ... , 

a ,<f> s series terminates if one or more of the numerator parameters is of the form q~ m for 
rn = 0, 1, 2, ... and q j* 0. When r = s + 1, the expression in (5) simplifies to 



s+\fts 



a \l > a s + 1 

b { , ... , b s 



; q ,x 



1 



n= 0 



(^1 i q)n ••• (*U+l > q)n 
(q ; q)n {b x ; q) n ... (b s ; c/)„ 



Note that in a basic hypcrgcometric scries £ c„jr, the ratio c n+ Jc n is a rational function 
of q n . [Ref. 12: p. 4] 

Using the q — shifted factorial (3) and (4), we can define basic extensions for many 
of the functions and formulas introduced in earlier chapters. We note that often there 
is more than one way to extend a result; examples will be given shortly. 

The q-gamma function is defined by 

(6) I%(.r) = — ' (1-?)'"*. 0<q< 1 . 

(<7 ; <7)«, 



Gosper showed that 

(7) lim. r („v) = r(x) 

<?- 1 

for 0 < q < 1 [Ref. 30: p. 109|. The structure of the gamma function extends as well. 
For instance, the formula 
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r,(.t+i) = 



\-q 



r,(.v) 



can be reduced using (6) to 



r(x+l) = xr(.r). 



With (6) we can define the q-beta function 



B q {x,y) 



r ? U) r q ( y) 
r q (x+y) 



which by (7) tends to Z?(.r, i y) as q -* 1". 

The q-binomial coefficient is defined for integers n and k by 




(<7 : <7 

(<7 i ?)* 0? I <?)„-* ’ 



where k = 0, For nonintcgral a and /?, we define 



'«! _ (/ +1 ; iL (?* /?+1 ; ?)oo 

_/?J? {q ; < 7 )^ (< 7 a+I ; < 7 )^ 

r,(a+l) 

" r,(/?+i)r ? («-/?+i) • 



The q-binomial theorem is then 



n 

{ab -,q) n = ^ j 



k = 0 



n 



b {a ; <74 (6 ; < 7 )„_* 



where n = 0, 1, 2, ... . [Ref. 12: p. 20| 

The next two expressions are basic extensions of the Chu-Vandermonde formula 
(Chapter III, Section F. 1) and are both known as the q- Chu-Vandermonde formula: 



2'Al 




c 



b 




(fib . ; </)„ 
(c ; < 7 ),! 
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2^1 



' —n , 

q , b 



c 




[Ref. 12: p. 11). These forms can be shown to be equivalent by reversing the order of 
summation. 

Jackson introduced the general form of the q-integral 



The q — integral defines the measure d q t which is a natural object for 



play an important role in generalizing the ordinary orthogonal polynomial classes. 

D. BASIC EXTENSIONS OF ORTHOGONAL POLYNOMIALS 

In this section, we present the q — analogue(s) of selected classes from earlier chap- 
ters. By using formulas such as those presented in the previous section together with 
methods based on those outlined in the preceding chapters, it is possible to derive the 
recurrence relations, difference equations, and Rodrigues' formulas as well as many 
other identities satisfied by these q —versions. As mentioned in the previous section, the 
q — extension of a function is not necessarily unique; however the last two classes listed 
are especially important. 

1. Continuous q — Herniite Polynomials 




where 




oo 



q — defined functions [Ref. 12: p. 19). All of the functions and formulas developed above 



a. Definition 




where x = cos 0. 
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b. Orthogonality Relation 



[ 



o H m { cos 0 1 q) ll n { cos G\q) \ (e W | 2 dQ = S m,n 



where .v = cos 0. [Ref. 12: p. 188] 

2. Discrete q — Hermite Polynomials 
a. Definition 



0 / 2 ] 



«„(*;?) = E 



(<7 ; q) n 



k= 0 if > 9 )/( (? > O^n-lk 



, *(fe-l) /t— 2 * 

(- 1 ) <7 '* 



A. Orthogonality Relation 



J* a m (x ; < 7 ) H n {x ; ?) d\jj{x) = q( 2 ) (9 ; <?)„ <5 mn 



where ^(x) is a step function with jumps 



l x 1 (* V ; </ L (‘i ■ ‘i ), 



(? 2 ;? 2 ) 



at the points x = ± q > , j — 0, 1, 2, ... . [Ref. 12: p. 193] 
3. q — Laguerre Polynomials 
a. Definition 
For a > — 1 , 



, W( . (?’ + ‘-w). 



i^i 



9 rt+2+l 

2 - 1- 1 > *7 > X( l 

q 



b. Continuous Orthogonality 



00 2 r 

L mi x '■ </) i 7) ( n ' V 7 l'\ -x 
0 (— (I— <7Jx , </) c 



r(«+i) r(-«) (V +1 ; <7),, 
r«(-a) 0 / ; q)n f /" 



m,n 
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c. Discrete Orthogonality 



E L £ > ( c< ? fc '* r) L n\ c( i k ; <z) 

k~~oo 



k(*+ 1 ) 

2 

(-c( !-</>/*; < 7)00 



(^ +l > </)« 

/ •. /! (> m,, 

(< 7 ; Rh R 



where 



(R •' <7)00 (~c(l~< 7 )< 7 a+1 ; <7)00 (-(1 -q)l{cq*) ; <7)00 
(?' * +I ; <7)00 ( _c ( 1 -R) ;< 7 )oo (--f( 1 -^w)oo 



[Ref. 12: pp. 194-195] 

4. Little <7 —Jacobi Polynomials 
Definition 



p n (x ; a, 6 ; < 7 ) = 2 0 , 



' — A 2 , 77+1 

R , 

aq 



\r, xq 



b. Orthogonality Relation 



OO 

E ; a. * ; ?) P„(<7* \a,b;q) 

x=0 



(bq ; ?) x 
(<7 ; r) x 



{aqf 



[q ; <?)„ {l-abq) ( bq ; <7),, (afy 2 ; <7)^ 
(<% ; q) n ( l-abq 2n+l ) (a? ; q) n {aq ; <7)^ 



{aqf 5 



m,n 



where 0 < q , a</ < 1 [Ref. 12: p. 166] 
5. Big q —Jacobi Polynomials 
a. Definition 



P n {x ; a,b,c\q) = 3 0 2 




, abq 



n + 1 



aq, cq 



x 

\q> <1 
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b. Orthogonality Relation 



ra< 
^ eg 



aq {x!a ; eft™ ( x l c y q)™ 

P m (x ; a,b,c-q ) P n (x ; a, b, c ; q) V - V 

C* ; q) oo ( bxlc ; <7)^ * 

,, {q\q) n ^-abq){bq\q) n {abqlc\q) n f'"') c 

= A* ; TTTT (— ac ) (7x2/ 

(a^ ; </)„ (l-afy ) («<7 ; (c<7 ; q) n 



where 



f a<? (Jf/fl ; <7)oo We ; ?)qq 

l f<? (■* ; <7)00 ( bx l c ; <7)00 ** 

gq_ ( i-< 7 ) (<7 ; ?)oo (fjf ; gU ( a ?/ c ; <7)00 ; <7)0, 

H ; ^)oo ( b q ; <7)00 ( c< 7 ; <7)00 ( ab ql c ; <7)00 



[Ref. 12: pp. 167-168] 

6. <7 — Krawtehouk Polynomials 

a. Definition 



K n (x ; a, N 4 , q) = 3 4> 2 



~ —n —l/i 

q , X, -a q 

-N A '» < 7 > 
q , o 




Z». Orthogonality Relation 



N 

X ’» a ’ # ; <?) Kn(q~ x ; w ; < 7 ) 



jc=0 



(<7 ;<7).t 

(<7 ; <7)x 



(-«/ 



1 _i x * -^ 1N ] (<7;<7)i,(i + a ')(-« '<7 ;</)/, 

— ( q a '<q)N a q ^ 2 ' , _i w< , -1 2«w \ 

(-a ;<7j„(l + a <7 ) [q \q) n 

w ( /V+l\— /? rt(rt+l) r 

x (-a? ) <7 <5m,/i 



[Ref. 12: p. 185] 

7. q —Hahn Polynomials 
a. Definition 



Qn( x » A/ , ^/) 3^2 



“ — / 77-fl — J C 

, abq , </ 



9 



— iV *9* 9 
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b . Orthogonality Relation 



N 

Yj Q m (- x ; jV ; </) <?„(* ; «. w ; q) 



x=0 



W ; <?)* ; ?)/v-.x , r x 

{q\q)Ar,q) N - x ^ 

(abq 1 ; < 7 )^ (<7 ;?)„(! -(%) (^<7 ; q) n (abq lV+2 ; ^)„ 



(<7 ; q\v {abq ; ^)„ (l -abq 2n+l ) (aq ; ?)„ (q A ; ?)„ 

(-^^( 2 ) Nn 5 mn , m, rt = 0, \, ... , N 



[Ref. 12: p. 165] 

8. Dual q —Hahn Polynomials 

a. Definition 

R n { q{x ) ; b,c,N ; q) = 3 <j> 2 

where ^(x) = q~ x + cq x ~ N 

b. Orthogonality Relation 



r — n —X X — <V 

q , q , cq 

-/V , ■* q> q 

q , beq 



N 

I 

x=0 



R m (h(x ) ; b, c, N ; < 7 ) /?„(At(x) ; 6 , c, iV ; q) 



(cq ,V ; c/)^ ( 1 -cq 2x ,V ) ( 6 c <7 ; q) x (q N ; < 7 )* 

— — _ iV , _ iV ' q ^ (—beq) 

(q,q)x{i-cq ){b q \q) x \ \q)x 
(i lc-,q)v (q ; q) n (bq ; q) n 



( bq ' q)s {beq ; q) n ( q A ; q) n 



{cq A )" b m n , m, n = 0, 1, , iV 



[Ref. 12: p. 166] 

9. ^ — Racah Polynomials 

a. Definition 



( 8 ) 



; a, b, c, d ; q) = 4 0 3 



— 7 ? / 72+1 —X / ,T +1 

q , abq , q , caq 



aq, bdq, cq 



> */> q 



where //(x) = q- x + cr/</' +l 
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b . Orthogonality Relation 

When aq , bdq , or cq = < 7 “* for a positive integer N, then 



N 

(9) 'Yj PM X ) 5 b > c > d 5 ?) AiG*M ;a,b,c,d; q) w(x ■ a, b, c,d \ q) = h n S m n 

x=0 

where 



w(x ;a,b,c,d; q) 



( 10 ) 



; gk W ; ; <?)jc (c? ; ?)x 

(<7 ; <7)* (« -1 ^ ; <i) x { b ~' cc i ; <7)* (<ty ; q) x ( abc if 



1 -cV ^ 1 

1 — cdq 



and 



(ID 



{q ; ?)„ (l-afrq) (fr<7 ; q) n ( acT^ q ; q) n ( abc ~ l q ; ?)„ («/<?)" 
(^<7 ; q) n ( \-abq 7n+X ) (a<7 ; q) n ( bdq ; ?)„ (cq- ; ^)„ 
H? 2 ; ?)oo ; <7)^ (a~ x d\ q)^ (fr -1 ; <7)00 

(« -1 c</<7 ; <7)^ (^ _1 c<7 ; q ) 00 (<fy ; <7)^ (a - ^ -1 ? -1 ; <7)00 



[Ref. 20: p. 1014]. When aq , bdq , or cq is equal to q~ N , the infinite products in (11) 
reduce to finite products. Hence the orthogonality relation (9) is valid for all q provided 
no zeros are introduced into denominator terms [Ref. 18: p. 4]. 

10. Askey-Wilson Polynomials 

a. Definition 

p n (x \a,b,c,d | q) = a~ n ( ab ; q) n (ac ; q) n {ad ; q) n 

(12) q , cibcdq , ae , ae 

X 403 , , !?1 <7 

a/;, ac, aa 

where x = cos 0 [Ref. 18: p. 3] 

b. Orthogonality Relation 
For —1 < a, b, c, d , q < 1 , 



1 f 1 

" 2 ^" J Pm{ x \a,b,c,d | < 7 ) t?„(x ; a, b, c, <-/ I < 7 ) 



tv(x ; a, b, c, d | q) dx 



1 -x 



b n b m,n 
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where 



Y\(l-2(2x 2 -l)< 1 k + q 2k ) 
w(x ;a,b,c,d\ q) = /,(*,«) h(x,b) h(x,c) h(x,d) 



with 



h(x,a) = J (l — 2axq k + a 2 q 2k ) 
k = o 



and 



( abcdq_ 2n _ ; q)*> ( abedg " 1 ; q) n (?" +1 ; <7) J (>V ; ?)J 
; <?)«, OV ; <7)oo ; <7)oo (W?" ; < 7 )^ {cdq n ; < 7 )^ ' 



[Ref. 18: pp. 11-14] 

We can establish a formal connection between these last two classes. If the 
q — Racah polynomials are written 



403 



" — rt , 1/ rt+1 - X , J, jc+1 

q , a b q , q , c d q 
a'q, b'd'q, c'q 



;<7> q 



and the following parameter changes are made: 



a = 



ac 

<7 





d' 



a 



b ‘ 



(or likewise using any permutation of { b, c, d) assuming at least one of these is non- 
zero), as well as the change of variable q- % = ae' 6 , then we obtain 
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—n , j n— 1 iO —iO 
q , abedq , ae , ae 

ab, ac, ad 



;?» q 



with 



q(x) = q- x + q x+l c’d’ 

10 —iO n 

= ac +ae = 2a cos (/ 
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That is, the q — Racah and Askey-Wilson polynomials are virtually the same, differing 
only in their parameters, normalization, and variable. 

Because the q — Racah and Askey-Wilson polynomials are essentially the 
same, the names are often used interchangeably in the literature. We have presented 
these polynomials as distinct classes in order to emphasize the continuous and discrete 
natures as expressed in the orthogonality relations. 

As stated at the beginning of Chapter IV, the "classical orthogonal 
polynomials" are defined to be those which are special or limiting cases of the Askey- 
Wilson (21) or q— Racah polynomials (18) [Ref. 31: p. 57]. We can now make this 
statement a bit more precise. Letting 

a= q , b = q H ,c = q y ,d=q 

in (9) and (10) and taking the limit as < 7 ~+l, the ordinary Racah polynomials 
R„( 2(.v) ; a, /?, y, 6) and their weight function as defined in Formulas (30) and (32), 
Chapter IV, Section C.3 are retrieved. A similar limiting process will recover any ordi- 
nary orthogonal polynomial class from its q —extension. Moreover, any < 7 —extension 
previously discussed is a special or limiting case of the q —Racah. For example, letting 
cq = q~ N and d = 0 produces the q— Hahn, etc. In this way, we see that the self-dual 
q— Racah polynomials encompass all the previous classes. They also satisfy three-term 
recurrence relations, second order difference equations, and Rodrigues' formulas with 
respect to < 7 —divided difference operators. The interested reader is referred to [Refs. 18, 
12 , 20 ] for details about these very rich classes. 

C. CONCLUDING REMARKS 

In Chapter V, we presented a few of the traditional applications in which orthogonal 
polynomials have played an important part. More recently, an enriching interplay has 
developed between the theory of orthogonal polynomials and other mathematical and 
mathematically-related areas. 

Efficient computational methods have been devised for determining the many useful 
quantities associated with orthogonal polynomials (such as their zeros, recurrence coef- 
ficients, etc.) [Ref. 32: pp. 181-216). Various classes of orthogonal polynomials have 
also played a role in digital signal processing [Ref. 32: pp. 1 15-133], quantum mechanics 
[Ref. 32: pp. 217-228], and birth/death processes [Ref. 32: pp. 229-255]. Advances in the 
field of combinatorics and graph theory have allowed new geometric interpretations of 
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orthogonal polynomial identities, some of which have very important consequences for 
"association schemes" and the designs of codes [Refs. 33, 16, 34, 32 : pp. 25-53, 35J. 

Physicists have introduced various versions of "diagrammatic methods": ways of 
understanding orthogonal polynomials through 3-j and 6-j symbols, their generaliza- 
tions, and accompanying identities by formally associating them with pictorial schemat- 
ics representing forces of physical systems that conserve angular momentum [Ref. 36]. 
Powerful new techniques involving "quantum groups" have been used to generate new 
identities for some classes [Ref. 32: pp. 257-292]. Further investigation into the 
electrostatics problem discussed in the text has led to the formation of the famed 
"Selberg beta integral" and its generalizations. It has yielded to analysis via the study 
of the "root systems" of Lie algebras, and has found applications ranging from statistical 
mechanics to computer algorithm complexity [Refs. 30: pp. 48-52, 32 : pp. 311-318, 37). 

Finally, research into the general structure of q —series has led to many surprising 
connections, and is intimately related to the many remarkable and powerful number- 
theoretic formulas discovered by S. Ramanujan, the famed Indian mathematical prodigy 
[Refs. 30: pp. 87-93, 38: pp. 55-66, 39]. 

This utility and promise of future applications provide ample motivation and justi- 
fication for continued study of the intrinsic structure of orthogonal polynomials. 
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